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20.  Abstract  (Continued) 


A  detailed  investigation  is  undertaken  of  the  effects  of  temperature  fluctuations  in  the  ocean 
on  the  phase  stability  of  the  signal  beam.  Both  large-scale  temperature  fluctuations  caused 
by  internal  waves  and  small-scale  fluctuations  caused  by  velocity  turbulence  and  thermal 
diffusion  are  considered.  Numerical  calculations  are  made  of  specific  cases  of  interest.  It  is 
concluded  that  the  stability  of  phase  is  either  strongly  affected  or  moderately  affected, 
depending  on  the  length  of  the  fiber  exposed  to  the  medium,  the  depth  in  the  ocean,  the  track 
length  of  towed  sensors,  the  time  of  exposure  in  moored  sensors,  and  the  characteristics  of 
built-in  instrumental  spectral  windows. 

Fart  II: 

The  effect  of  extraneous  (environmental)  mechanical  stress  on  the  state  of  polarization  of 
light  signals  in  fiber-optic  sensors  of  acoustic  pressure  is  investigated.  After  review  of  the 
relevant  theory  of  the  propagation  of  polarized  light  in  a  dielectric  waveguide  and  of  photo¬ 
elastic  coupling,  more  specific  cases  of  the  birefringent  effects  of  mechanical  stress  are 
theoretically  analyzed  and  numerically  calculated.  It  is  concluded  that  these  effects  are 
negligible  in  the  cases  treated  except  for  fibers  subject  to  equal  and  opposite  forces  along  a 
diameter  or  subject  to  random  pressure  fluctuations.  Estimates  of  differentia)  changes  in 
the  index  of  refraction  in  silica  fibers  due  to  (elastic)  stress-induced  birefringence  are  furnished. 

;n 


-  "1 

i 

X 


1 


U  73fi 


security  classification  of  this  PAGE'**r>*n  n»t»  > 


j 


CONTENTS 


SUMMARY  . 

INTRODUCTION  . . . 

CHARACTERIZATION  OF  THE  STATE  OF  POLARIZATION  OF  LIGHT . 

TYPES  OF  BIREFRINGENT  MEDIA  . 

EVOLUTION  OF  THE  STATE  OF  POLARIZATION 
WITH  DISTANCE  OF  PENETRATION  . 

EFFECTS  OF  DIELECTRIC  POLARIZATION  ON  THE  PROPAGATION  OF  LIGHT 
IN  A  WAVEGUIDE  . 

Field  Equations . 

Mode  Coupling  in  a  Dielectric  Waveguide  . 

Dielectric  Polarization  by  Torsional  Stress  . 

Mode  Coupling  Coefficients  for  Pure  Torsion . 

Evolution  of  Electric  Vector  Polarization 

in  the  Case  of  Pure  Torsion  of  a  Fiber  . . 

Periodic  Dielectric  Polarization  . 

ANALYSIS  BY  PHOTOELASTICITY  OF  ELECTRIC  FIELD  POLARIZATION  OF  SINGLE 
MODE  FIBERS  UNDER  MECHANICAL  STRESS  . 

Stress  Analysis . 

Propagation  of  Polarized  Light  in  Doubly  Refracting  Materials  . 

Polarization  Effects  in  Stressed  Elastic  Bodies . 

Summary  . 

APPLICATIONS  . . 

NUMERICAL  CALCULATIONS  . 

Stress-Optic  Constant  C . 

Numerical  Calculations  of  Fiber  on  a  Mandrel  (Case  6)  . 

Bare  Core  on  a  Mandrel  (Case  6)  . 

Parameter  Variation  for  Calculations  of  Aw  of  a 

Bare  Core  on  a  Mandrel  (Case  6)  . 

Fiber  Squeezed  by  Diametral  Press  (Case  8)  . 

Average  A//  Along  a  Horizontal  Diameter  (Case  8)  . 

Pure  Torsion  (Case  9)  . . 

Fiber  Subject  to  Random  Stresses  (Case  8)  . 

CONCLUSIONS  . 

REFERENCES  . 


42 


SUMMARY 


Polarized  laser  light  conducted  in  optical  fibers  is  modulated  by  time-varying  and  space-varying 
mechanical  stress.  Stress  effects  on  stales  of  polarization  are  particularly  important  in  applications  of 
optical  fibers  to  sensing  of  acoustic  and  mechanical  fields.  Quantitative  estimates  are  essential  in  these 
cases.  The  procedure  for  making  these  estimates  followed  these  steps:  first,  to  lay  a  firm  groundwork 
for  making  numerical  calculations  a  detailed  review  of  the  theory  of  coupling  between  electric  fields  and 
mechanical  stress  was  undertaken.  An  important  distinction  was  made  between  single-mode  fibers  car¬ 
rying  one  mode  and  the  same  fibers  carrying  two  degenerate  modes.  The  results  of  this  review  are 
presented  in  a  group  of  equations  which  are  applied  to  the  case  of  pure  torsion.  The  resultant  calcula¬ 
tions  agree  with  the  work  of  Ulrich  and  Simon.  Upon  completion  of  this  review  a  number  of  specific 
cases  were  calculated  of  the  change  in  phase  and  state  of  polarization  of  polarized  light  in  optical  fibers 
disturbed  by  mechanical  stress;  these  cases  were 

Case  1.  Applied  stress  on  the  fiber  is  axisymmetric  and  everywhere  uniform  along  the  optical 
path 

Case  2.  Applied  stress  on  the  fiber  is  axisymmetric  but  varies  along  the  optic  path 

Case  3.  Tensile  or  compressive  stresses  are  applied  to  the  ends  of  the  fiber 

Case  4.  Oblique  incidence  of  optic  ray  and  pure  bending 

Case  5.  Fiber  subject  to  both  bending  and  shear 

Case  6.  Fiber  wrapped  on  a  mandrel 

Case  7.  Fiber  subject  to  random  stress 

Case  8.  Fiber  squeezed  by  diametral  forces 

Case  9.  Pure  torsion 

It  is  concluded  that  the  birefringent  effects  are  negligible  in  most  applications  because  there  are  no 
components  of  stress  transverse  to  the  optic  ray.  Two  important  cases  do  show  birefringent  effects, 
fiber  squeezed  by  diametral  forces  and  fiber  subjected  to  random  stresses.  Numerical  calculations 
predict  that  for  fibers  of  fused  silica,  at  the  limit  of  allowable  stress, 

Aflmax  5=5  “8.5  x  10-4  for  the  case  of  applied  diametral  forces 


and 

=  1  rad  if  the  rms  stress  fluctuation  is  23.4  kPa  (3.4  psi). 

The  case  of  a  fiber  80  /xm  in  diameter  with  a  core  4.5  jwm  in  diameter  wound  on  a  mandrel  1  cm  in 
diameter  has  also  been  numerically  estimated  for  fused  silica  at  the  maximum  allowable  stress.  Here 
the  birefringent  effect  is  negligible  if  the  optic  ray  is  parallel  to  the  axis  of  the  fiber.  However,  if  the 
optic  ray  is  incident  at  4°  to  the  optic  fiber  axis,  then  it  is  estimated  that 

A »MI  =  -2.5  x  10  \ 
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OCEAN  ENVIRONMENTAL  INFLUENCES  OF  TEMPERATURE 
AND  MECHANICAL  STRESS  ON  BARE  FIBER  OPTIC 
SENSORS  OF  ACOUSTIC  PRESSURES 
PART  II 


INTRODUCTION 

The  effect  of  an  optical  waveguide  fiber  under  various  conditions  of  temperature,  pressure,  etc. 
upon  the  state  of  polarization  of  incident  light  has  been  investigated  by  several  authors  [1-3].  Recently, 
these  studies  have  been  made  in  connection  with  the  magnetooptic  (Faraday)  effect  [4,5l.  Although 
varied  in  details,  most  papers  present  the  typical  laboratory  experimental  setup  of  Fig.  1:  A  He-Ne  laser 
(1)  generates  linearly  polarized  light  with  a  fixed  plane  of  polarization.  This  passes  through  a  quarter- 
wave  plate  (2)  and  emerges  circularly  polarized.  The  light  then  enters  a  polarizer  (a  Glan-Thompson 
prism)  (3),  which  is  mounted  on  a  rotatable  fixture,  and  emerges  linearly  polarized  at  any  desired  polar¬ 
izing  angle.  This  beam  is  then  focused  by  a  microscope  objective  (4)  onto  a  fiber  (5).  A  watertight 
compartment  (6)  allows  measurement  at  various  fiber  temperatures.  The  emerging  light  is  collimated 
by  a  lens  (7)  and  then  analyzed  by  a  rotatable  Wollaston  prism  (8),  which  splits  the  incident  light  into 
slow  and  fast  components  by  birefringence  (9).  These  components  are  detected  by  photodetectors  (10) 
which  deliver  intensities  /j  and  I2 


BATH  FLUID 

Fig.  1  —  Fxperimentul  stud>  of  polarization 


During  the  experiment  the  polarizer  (2)  is  set  to  an  arbitrary  angle  (input  polarization  angle) 
and  the  analyzer  (8)  is  rotated  to  an  output  angle  ^  such  that  one  of  the  emerging  beams  has  a  max¬ 
imum  intensity  /max  and  the  other  has  a  minimum  intensity  /min.  The  degree  of  polarization  P  is  then 
determined  by  analog  electronic  circuits, 

p  _  Anin 

^max  4"  ^min 

Thus  polarization  effects  of  a  fiber  of  length  L  are  characterized  by  measurement  of  <t>,  0,  and  P  as  a 
function  of  environment  (temperature,  pressure,  etc.).  A  birefringent  model  of  the  fiber  then  relates  L 
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to  the  phase  difference  8  between  fast  and  slow  modes.  In  this  way  8  is  deduced  by  measurement.  The 
magnetooptic  effect  on  polarization  is  studied  by  introducing  a  magnetic  field  with  a  coil  (11). 

In  applications  of  optical  waveguides  to  such  tasks  as  communications,  measurement,  and  control 
of  information  and  processes  the  concept  of  depolarization  is  crucial.  This  refers  to  reduction  in  P  and 
change  (possible  randomization)  in  0  and  8  due  to  temperature  and  mechanical  stress.  It  also  refers  to 
the  loss  of  coherence  between  fast  and  slow  waves  in  the  fiber  caused  by  the  same  external  factors.  A 
summary  of  facts  and  theories  of  partially  polarized  light  is  found  in  standard  treatises  [6,  p.  552], 

This  report  will  concentrate  on  the  theory  of  the  effect  of  mechanical  stress  or  strain  on  the  state 
of  polarization  of  a  light  beam  in  an  optical  waveguide  fiber.  An  experimental  setup  to  study  this  effect 
could  be  the  same  as  in  Fig.  I  except  for  replacement  of  the  temperature-controlled  bath  by  various 
strain-induced  deformations  of  the  fiber. 


CHARACTERIZATION  OF  THE  STATE  OF  POLARIZATION  OF  LIGHT 


A  brief  review  of  standard  theory  will  help  ciarify  more  advanced  concepts  in  this  report.  The 
propagation  of  polarized  light  E  *=  (Ex,  Ev,  Ez)e,Uit  in  a  fiber  can  be  visualized  as  the  trace  of  the  tip  of 
the  vibrating  electric  vector  transverse  to  the  direction  of  propagation.  In  general  this  trace  is  a  helix 
whose  radius  varies  with  distance  of  propagation.  The  projection  of  an  adjacent  portion  of  the  helix  on 
any  transverse  plane  is  an  ellipse  whose  semimajor  axis  is  a  units  long,  oriented  at  angle  $  relative  to 
the  x  axis,  and  whose  semiminor  axis  is  b  units  long.  The  helix  is  called  right  handed  (or  positive)  if 
an  observer  looking  face  on  at  the  source  of  light  sees  the  electric  vector  rotating  clockwise;  it  is  left 
handed  (or  negative)  if  the  vector  is  rotating  anticlockwise.  A  reversal  in  direction  of  the  ellipse  means 
a  switch  from  positive  to  negative,  that  is,  the  azimuth  angle  \  changes  by  n  radians.  To  any  ellipse 
with  its  major  axis  at  azimuth  a  there  corresponds  an  orthogonal  ellipse  whose  major  axis  is  at  90°  +  a. 
Since  the  ellipse  can  be  represented  in  parametric  (or  component)  form. 


and 


or 


Ex  =  a  cos  (r  +  8  () 
£v  =  b  cos  (r  +  S2), 


Ev  b  —  /  (fi  |  —  6  o 

%  =  ~a  e 


its  orientation  (angle  <//)  depends  on  8j  —  82.  When  the  phase  8  changes  with  distance  of  travel  z,  its 
magnitude  at  wavenumber  k  is 

8  =  kz  =  kQnz  =*  nz, 

where  n  is  the  index  of  refraction  and  X0  is  the  wavelength  of  the  light.  Thus, 

8  }  —  5  2  —  ~  z(n,  -  n2) 
k  o 

In  an  anisotropic  (birefringent)  fiber  n\  ^  n2.  This  causes  the  ellipse  to  change  orientation  as  it  pro¬ 
pagates. 

A  positive  (clockwise)  ellipse  at  fixed  azimuth  \  ranges  in  appearance  from  a  line  (6  —  0)  to  an 
ellipse  (b  ^  a)  to  a  circle  (b  =  a).  A  negative  ellipse  exhibits  the  same  series  of  shapes,  except  with 
reversed  rotation. 
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The  variation  of  the  components  Ex,  Ev,  Ez  propagating  in  an  optical  fiber  along  coordinate  z  can 
be  studied  by  the  matrix  algebra  of  the  Jones  calculus  [7].  An  equally  illuminating  procedure  is  to 
model  the  state  of  polarization  of  the  light  as  a  point  on  the  surface  of  a  unit  (Poincare)  sphere  .V  [8] . 
On  this  sphere  an  arbitrary  point  represents  a  state  of  elliptic  polarization  which  is  traced  out  in  time  by 
the  electric  vector  at  a  particular  cross  section  z  —  const,  of  the  fiber.  All  latitudes  above  the  equator 
represent  positive  rotation;  those  below'  represent  negative  rotation.  Since  the  relative  shape  of  the 
ellipse  can  be  given  by  the  ratio  b/a  of  its  axes,  and  its  orientation  by  an  angle  0  of  the  major  axis  rela¬ 
tive  to  some  reference  line,  the  point  on  S  can  be  located  by  two  angles:  an  azimuth  angle  2t //,  and  an 
elevation  angle  2\  (Fig.  2).  The  angle  appears  double  because  an  ellipse  completes  a  cycle  of  orienta¬ 
tion  in  a  plane  cross  section  in  v  radians  (180°),  whereas  a  point  P  completes  the  same  cycle  on  S  in 
271  radians  (360°).  Also,  since  \  varies  from  0  to  45°  as  b/a  varies  from  0  to  1,  it  is  seen  that 

tan  \  =  b/a. 

Two  special  cases  of  the  elliptical  slate  of  polarization  appear  on  S. 

Ellipse  degenerates  to  a  line  (b  =  0) 

Here  \  —  0,  hence,  all  polarization  states  P  are  found  on  the  equator  (Fig.  2).  At  point  H  the 
angle  2ii>  —  0;  the  light  is  then  linearly  polarized  in  the  direction  of  the  x  axis.  At  point  V  the  angle 
2\b  -  7r,  or  \b  =  90°;  the  light  is  then  linearly  polarized  in  the  direction  of  the  y  axis.  At  C  the  plane  of 
polarization  is  45°  relative  to  fhe  a  axis,  and  at  D  it  is  135°. 

Ellipse  degenerates  to  a  circle  (b/a  =  1 ) 

Here  tan  x  =  1  (i.e.,  2\  =  90°).  There  are  only  two  states  of  polarization,  L  and  R  (Fig.  2). 
Point  L  (north  pole)  represents  circularly  polarized  light  rotating  counterclockwise  (when  looking  from 
L  toward  the  origin),  while  point  R  (south  pole)  represents  circularly  polarized  light  rotating  clockwise 
(when  looking  from  R  toward  the  origin). 


Rig  2  —  Angles  and  2\  on  a  Poincare  sphere 


R 


When  a  light  beam  of  specified  polarization  (P \  on  S)  enters  a  birefringent  medium,  its  state  Pi  at 
increment  dz  of  penetration  is  a  point  on  the  trajectory  dC(z)/dz  of  the  unit  radius  vector  C(r)  on  S 
whose  initial  terminus  is  P{.  This  change  from  Pt  to  P2  is  caused  by  the  phase  advance  8  with  distance 
of  penetration  of  the  fast  component  of  the  light  vector  over  the  slow  component,  which  thus  alters  the 
initial  time  phase  of  the  light  vector  components  along  the  major  and  minor  axes  of  the  ellipse  at  P]. 
Since  any  incremental  trajectory  on  the  surface  of  a  sphere  corresponds  to  a  rotation  about  some  diame¬ 
ter,  one  can  find  state  P2  by  rotating  the  entire  sphere  S  through  an  angle  8  about  a  uniquely  chosen 
axis  of  rotation  AAU.  where  A  and  Aa  are  opposite  points  at  the  ends  of  a  single  diameter,  representing 
the  two  unique  polarization  stales  (or  eigenstates)  of  the  medium  relative  to  some  reference  line. 
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For  example,  assume  that  a  linear  birefringent  medium  is  so  oriented  that  it  resolves  an  incident 
linearly  polarized  light  into  two  linearly  polarized  orthogonal  components,  one  along  X  and  the  other 
along  K  Directions  X  and  Y  are  taken  to  be  the  privileged  axes  along  which  the  states  of  polarization 
do  not  change  with  penetration.  The  component  along  X  is  at  state  //  on  S ,  while  the  component  along 
Y  is  at  state  V.  The  axis  HV  is  therefore  the  axis  of  rotation  for  any  state  P\  moving  along  trajectory 
C(z)  to  state  P2  To  illustrate,  let  P\  be  a  linearly  polarized  beam  incident  on  the  birefringent  medium 
at  00  —  20°  relative  to  X.  Its  point  on  S  is  on  the  equator  at  azimuth  angle  40°  relative  to  HV.  It  is 
resolved  by  the  medium  along  the  privileged  axes  X and  F,  where  the  resultant  components  travel  at 
different  speeds.  After  the  beam  penetrates  a  distance  AZ,  the  phase  along  the  fast  axis  (say  //  ) 
advances  (say)  5  =  15°  ahead  of  the  phase  along  the  slow  axis  (say  V).  The  state  P2  (A Z)  on  .S  is  then 
obtained  by  rotating  the  entire  sphere  about  HV  counterclockwise  (looking  from  the  fast  axis  //  toward 
the  origin)  through  an  angle  of  15°  (Fig.  3).  Since  the  elevation  angle  2\  of  P2  is  different  from  zero, 
it  is  seen  that  state  P2  represents  elliptical  polarization. 


c'hn  0 

COS^o 

£> 

o 

sinP0 

The  angle  0  of  the  major  axis  satisfies  the  equation 

tan  20  =  tan20ocos8  =  tan  40  °cosl  5  0=s  0.81, 

or  20  *=  39.02°  and  0  —  19.51°.  The  relations  between  £max,  £mm,  £v .  £x,  £(),  0,  andfl0  are  shown  in 
Fig.  4.  The  new  coordinates  of  polarization  state  P2  are  azimuth  angles  20  2  and  elevation  angles  2x:. 

In  sum,  in  this  example  the  incident  wave,  linearly  polarized  at  f>()  =  20°,  becomes  elliptical Iv 
polarized  at  19°.  The  magnitudes  of  £max,  £min,  an^  degree  of  polarization  will  be  discussed  later. 


L 


t  ig.  3  —  t  volution  of  suite  nl  pol.m/.ition  from  P,  to  P-. 
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pig.  4  -  Definitions  «>t  parameters  in  elliptical  polarization 


TYPES  OF  BIREFRINGENT  MEDIA 

The  axis  of  rotation  AAa  on  the  Poincare  sphere  is  determined  by  the  optical  properties  of  the 
medium.  Three  types  of  media  can  be  described. 

A  medium  which  has  two  eigenstates  of  linear  polarization  is  called  a  linear  birefringent  medium. 
When  a  beam  with  elliptical  polarization  is  incident  on  this  type,  it  is  resolved  into  two  orthogonal 
linearly  polarized  beams  (referred  to  as  the  fast  beam  and  the  slow  beam).  In  this  medium  all  axes  of 
rotation  are  diameters  confined  to  the  equatorial  plane  of  the  S  sphere.  The  unique  diameter  (axis  of 
rotation)  for  a  given  state  P]  is  fixed  by  the  orientation  of  the  privileged  birefringent  axes  relative  to 
the  Tand  Y  axes  of  the  (transverse)  coordinate  system.  In  the  above  example  the  unique  diameter  w'as 
H V.  The  direction  of  rotation  is  always  taken  to  be  counterclockwise  around  the  faster  state  when  look¬ 
ing  toward  the  origin. 

A  medium  which  has  two  eigenstates  of  circular  polarization  is  called  a  circular  birefringent 
medium.  When  a  beam  with  elliptic  polarization  is  incident  on  this  type  of  medium,  it  is  resolved  into 
two  circularly  polarized  beams  of  opposite  rotation.  The  unique  diameter  of  rotation  is  the  axis  LR  on 
the  S’ sphere.  The  trajectory  of  new  states  ip2(AZ)  lies  on  latitude  circles,  and  the  direction  of  rotation 
is  counterclockwise  about  the  faster  state,  as  in  the  previous  case.  When  the  incident-beam  polarization 
is  linear,  this  medium  resolves  the  beam  into  two  oppositely  rotating,  equal-amplitude,  circularly  polar¬ 
ized  beams.  The  net  effect  after  the  beam  penetrates  a  distance  AZ  is  rotation  of  the  incident  plane  of 
polarization  to  a  new  plane.  This  medium  exhibits  optical  activity. 

A  medium  which  has  two  eigenstates  of  elliptical  polarization  is  called  an  elliptically  birefringent 
medium.  When  a  beam  with  elliptic  polarization  is  incident  on  this  type  of  medium,  it  is  resolved  into 
two  elliptically  polarized  beams  with  opposite  rotation.  The  unique  axis  of  rotation  for  a  given  state  P\ 
is  an  arbitrary  diameter  of  the  S  sphere.  This  medium  exhibits  a  combination  of  birefringence  and  opti¬ 
cal  activity.  The  incremental  rotation  A'  due  to  an  incremental  birefringent  phase  difference  8'  and  an 
incremental  optical  activity  p'  is 

A'  =  Vfi'2  +  (2p')2. 

The  unique  axis  of  rotation  for  a  state  P\  is  a  diameter  lying  in  the  plane  of  the  corresponding  linearly 
birefringent  medium  axis,  but  at  an  elevation  such  that 
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arctan 


V 

8' 


The  direction  of  rotation  is  counterclockwise  about  the  faster  state. 


EVOLUTION  OF  THE  STATE  OF  POLARIZATION  WITH  DISTANCE  OF  PENETRATION 

Ulrich  and  Simon  [9]  have  studied  the  effect  of  twist  on  the  state  of  polarization  of  single-mode 
fibers.  Their  approach  is  reviewed  here. 

Let  viz)  be  a  vector  in  xvz  space  whose  magnitude  teal  is  chosen  to  represent  the  angle  of  rota- 
lion  per  unit  length  of  penetration  and  whose  associated  spherical  angles  2*^  and  2£w  locate  the  end 
point  of  a  diameter  of  an  5  sphere  which  forms  the  axis  of  rotation.  Let  C (z)  be  a  unit  vector  of  this 
sphere  extending  from  the  origin  of  S  to  the  state-of-polarization  point  P\  at  distance  Z  along  the  fiber 
Then  the  new  state  of  polarization  a  short  distance  beyond  Z  is  given  by 


Assume  that  state  P\  represents  elliptical  polarization  and  assume  that  the  medium  rotates  P \  through 
an  angle  piz)  due  to  birefringence  and  an  angle  a(z)  due  to  optical  activity.  Then  the  total  rotation 
vector  for  small  increments  in  C  is 


to  (z)  =  a  (z)  +  0  (z). 


First,  let  <o,  a,  and  p  be  independent  of  z.  There  are  then  two  points  (±  co/co),  at  the  ends  of  a  fixed 
diameter  on  a  Poincare  sphere,  which  are  the  privileged  directions  of  propagation  of  the  two  eigenstates 
of  elliptic  polarization.  In  the  cross  section  (z  =  const.)  of  an  optical  fiber  the  orientation  of  these 
ellipses  is  given  by  p  (a  constant).  The  length  ratio  of  the  axes  is  given  by  the  elevation  angle  2\. 
where  fas  before) 


ja  1 

i,n*-  w 


Next,  let  a  and  p  depend  on  travel  distance  z  This  is  the  case  of  a  twisted  birefringent  fiber.  As  the 
polarized  light  penetrates  the  fiber,  two  events  occur  simultaneously:  (a)  the  privileged  axes  of 
birefringence  rotate  and  (b)  the  phase  of  the  fast  wave  relative  to  the  slow'  wave  increases.  Both  effects 
can  be  modeled.  To  model  event  (b),  Ulrich  and  Simon  [9]  insert  a  local  Cartesian  coordinate  system 
R°  at  each  z  =  z°,  and  orient  .v°  and  to  be  parallel  to  the  local  fast  and  slow  axes  respectively.  Polar¬ 
ization  changes  in  this  local  system  are  represented  on  the  5°  sphere  by  the  unit  vector  Cjj  rotating 
about  the  axis  of  rotation  H°  fixed  in  S°.  To  determine  the  magnitude  and  direction  of  O0,  one  must 
be  given  the  birefringence  rotation  /3°,  the  optical  activity  rotation  a0,  and  the  angle  of  twist  rotation  r°, 
all  in  R°.  Then 


n°  =  0°  +  a°-  2 r° 


The  factor  2  appears  here  because  azimuths  on  a  Poincare  sphere  are  twice  those  of  real  space  The 
minus  sign  is  a  convention  that  the  twist  is  clockwise  as  z  advances  while  the  rotations  p(1  and  a"  are 
counterclockwise.  In  the  local  system  diameter  ft n  defines  two  eigenstates  of  elliptical  polarization 
±n°/fl  which  propagate  in  R()  unchanged. 

To  model  event  (a),  the  entire  local  5°  is  rotated  about  LR  at  a  rate  2 r  relative  to  fixed  laboratory 
coordinates.  Thus  the  total  trajectory  C(z)  of  the  polarization  state  is  the  vector  sum  of  a  trajectory 
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along  the  rim  of  the  base  of  a  cone  on  .S’0  (apex  at  the  origin)  plus  the  trajectory  imparted  by  the  gen¬ 
erator  of  the  cone  as  its  tip  moves  in  a  latitude  circle  around  RL  of  the  laboratory-fixed  sphere  .V  The 
combination  makes  C(r)  trace  out  a  cycloid  on  S 

This  analysis  closely  parallels  that  of  the  general  motion  of  a  spinning  top  in  which  ft"  is  the  ana¬ 
log  of  the  angular  velocity  of  the  top  in  a  local  coordinate  system,  and  r  is  the  angular  velocity  of  the 
local  system  relative  to  the  fixed  laboratory  system.  The  top  undergoes  precession,  which  is  the  analog 
of  the  birefringent  effect  alone,  and  nutation,  which  is  the  analog  of  the  combination  of  optical  activity 
and  birefringence. 

The  trajectories  described  above  depend  on  specifically  located  axes  of  rotation.  The  determina¬ 
tion  of  these  axes  is  sufficiently  complicated  in  the  general  case  to  require  a  complete  review  of  the 
theory  of  the  effects  of  perturbations  of  the  dielectric  polarization  on  the  propagation  of  light  in  a 
dielectric  waveguide.  This  is  taken  up  next. 


EFFECTS  OF  DIELECTRIC  POLARIZATION  ON  THE  PROPAGATION 
OF  LIGHT  IN  A  WAVEGUIDE 

Field  Equations 

(1) 

(2) 

(3) 

(4) 


Maxwell's  equations  in  the  SI  units  form  the  starting  point  for  this  analysis: 

dD 


V  x  H  -  J  + 


dr 


and 


VxE+y=0. 

dt 

V  •  D  =  p, 


S7  ■  R  =  0 


where  the  units  are  H,  A  m  l;  J,  A  m  2;  D,  C  m  2;  E,  Vm  1 ;  C*m  3;  and  B,  V-s-rrf 

To  this  set  one  adds  the  constitutive  equations  which  relate  B  and  H  and  which  relate  D  and  E: 

B  =  MoH  +  /xoM-p'H,  fx'  =*  /x0(l  4-  MH  ')  (5) 

and 

D  =  e0ETP,  (6) 

where  the  units  are p ,  N'S2/C2;e,  C2/N  m2:  M,  A  m  l;  and  P,  Cm  2.  The  vectors  P  and  M  are  the 
dielectric  polarization  and  magnetization  respectively. 

Assume  J  =  0.  Elimination  of  B  from  (1)  and  (2)  leads  to  an  equation  in  E, 

V:E  —  V(V  •  E) -**'£„  -p’  -^7-  =  0.  (7) 

dr  dr 

In  this  analysis  we  neglect  magnetization  and  set  g '  -  g ()  ®  0(l)Ns2C  2  Also,  in  the  absence  of  true 
charge  p  the  displacement  D  is  a  constant.  Thus, 

V  D  =  V  -(€()E  +  P)  =  0  (8) 
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or 

V  *  E  =  — €0'!  (Vt0  •  E)  —  €0  1  V  P 

We  apply  Eq.  (7)  to  a  fused  silica  waveguide  of  very  small  diameter.  In  the  normal  stale  of  this 
material  we  assume  P  to  be  negligible.  In  the  perturbed  state  caused  by  mechanical,  thermal,  or  optical 
stresses  the  vector  P  takes  on  significant  magnitudes.  A  simple  model  which  is  useful  for  our  purposes 
has  P  proportional  to  E, 

P  =  UE,  (9) 


in  which  the  SI  units  of  the  elements  of  dyadic  U  are  those  of  the  dielectric  constant  e,  namely, 
C2/m2*N.  Applying  (8)  and  (9)  in  the  steady  state  (frequency  /  =  o)/2it)  to  (7),  one  obtains 

V2E  +  V  fco"1  V € y  •  E  +  € o  1  V  *  UE]  +  *2E  +€0-l*2(UE)  -  0,  (10) 

where 

/x'aj2€()=^^.0€0aj"  — *  k2. 

For  high  purity  fibers  we  can  readily  assume  the  €0  is  negligible. 

To  apply  Eq.  (10)  we  must  deal  with  components  of  vectors.  A  useful  decomposition  is,  in 
tangential  (superscript  /)  and  longitudinal  (superscript  z)  components, 

E  *  E'  +  k  Ez, 

V-<vV+£. 

V'  = 

ox  dy 
and 

P=  P'+  k  P:  —  (UE)'+  k( UEK 


In  this  report  interest  will  be  focused  mainly  on  the  tangential  vector  E'.  Its  propagation  is  governed  by 
the  reduced  equation, 

V2E'  +  V'(€0_1  v  •  UE)  +  A:2E'  +  €0-‘A:2(UE)'=  0.  (11) 

We  seek  a  solution  of  (11)  as  a  superposition  of  waveguide  modes  in  a  circular  fiber, 

E '(x,y,z)  *=  £  a„,(z)e'k m‘  E'm  (x,^)  (12a) 

m 

and 

U E (x,y,z)  =  J  o/M(2)(AZUE^,+  k(e^ZyJE h)2  (12b) 

m 

Now 

v2e  =  £  flm(z)Ar[(v')2 e;  -  kl e;i  +  £  k;  +  2 ikmam),  (13) 

m  m 

in  which  the  prime  in  am  means  9/9..  When  (13)  is  inserted  into  (11)  a  number  of  terms  of  the  latter 
cancel  because  they  satisfy  the  steady-state  wave  equation  of  waveguide  modes,  in  which  the  expansion 
coefficient  am(z)  —  am  =  constant. 
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(V')2E,',  +  (A:2  -  k*)E'm  =  0. 


(14) 


We  therefore  consider  only  the  second  sum  in  Eq.  (13).  In  it,  all  terms  in  ^"(z)  are  negligible  for  z  > 
\,  \  =  2-rr/k.  Thus  the  propagation  of  E,',  is  governed  by  the  equation 


ik„,: 


+  M 


m 

-,H  a2 


2  ikm  a’m  (z)UE£  +  jE>,)  +  am^'  A2[/(UEJ*  +  y(UEm)‘] 


(15) 


dx2(UEJV+ 


(UEm)' 


+  7 


— — (uF  )x  4-  — _ (LJ  F  )■' 

bydx  +  By2  m 


+  1 


Bxdz  m  JBydz  m) 


=  0. 


In  this  infinite  sum  one  can  find  a'm  in  terms  of  am  by  use  of  the  orthogonality  of  modes  in  the 
waveguide.  The  relation  a'Jam  determines  the  coupling  of  modes  caused  by  the  dielectric  polarization. 
Its  calculation  is  described  next. 


Mode  Coupling  in  a  Dielectric  Waveguide 

The  time-averaged  flux  of  energy  in  a  waveguide  is 

W  =  <f(E  x  H  *)  kdA  =  jj(E'  x  H  "Yd A  (16) 

For  each  mode, 

I'YdA.  (17) 


We  apply  this  formula  to  a  circular  fiber  whose  diameter  is  so  small  as  to  permit  the  propagation  of  a 
single  mode  only.  Since  this  mode  can  be  arbitrarily  oriented  in  a  cross-section  of  the  fiber,  it  can  be 
split  into  two  orthogonal  modes  propagating  with  the  same  magnitude  of  wavenumber.  These  are 
designated  as  degenerate  modes. 


Application  of  (17)  to  a  case  of  degenerate  modes  for  a  single-mode  waveguide  requires  special 
caution.  Assume  a  pair  of  degenerate  modes  with  imaginary  z  components  (this  example  is  used  by 
Ulrich  and  Simon  [9]) 


Ei  —  — 


M  H2 


Vm o 


* 0 


-  {ff.  E\,  £f }  - 


Jo(r),0,—  COS 
k\ 


(18a) 


and 


[El  Ei  E\ )  = 


—  sin 
k  i 


(18b) 


where  d>  is  a  polar  angle  in  the  cross  sectional  area  A ,  J0(r)  =  dJJdr,  k\  is  the  wavenumber  in  the  first 
(degenerate)  mode,  and  the  components  are  nondimensional.  If  (18)  is  inserted  into  (17)  and  the 
latter  is  applied  without  modification,  one  obtains 


wm  =  E'm  +  jE>„)  X  ( ///*  +  j  H-,)]dA 

-  hmn§(E*m  Hr  -  Et„  H't  *)dA. 


(19) 
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Since  the  goal  of  this  analysis  is  to  find  the  ratio  a'Jam ,  we  are  required  to  reduce  the  sum  on  a‘m(z)  to 
one  term.  We  therefore  multiply  the  /  component  of  (15)  by  //,j  *  and  the  j  component  by  //,;*.  When 
m  ^  n,  Eq.  (19)  is  zero  as  required.  However,  if  m  =  1  and  n  =  2,  it  is  seen  from  (18a)  that  H,\*  is 
zero.  Hence  all  terms  in  the  /  component  will  be  zero.  We  will  therefore  not  be  able  to  use  this  sum  to 
find  a  relation  for  a'm  in  terms  of  am.  On  the  other  hand,  H\  is  not  zero,  thus  the  j  component  can  be 
used.  To  avoid  this  difficulty,  we  use  a  second  orthogonality. 


Wmhmn=  X  ELVdA 

=  $WEH-  HiEr)dA. 


Again,  when  m  =  1  and  «?  —  2,  one  multiplies  the  /  component  of  (15)  by  H\.  Since  this  is  not  zero, 
the  /  component  will  yield  a  relation  between  a„  and  am .  In  contrast,  multiplication  of  the  i  component 
of  (15)  by  H\  causes  all  terms  to  vanish.  These  thoughts  lead  to  the  following  summary:  to  apply 
orthogonality  of  modes  to  (15),  we  multiply  the  /  component  by  and  the  j  component  by  //„\  and 
then  integrate  over  the  cross-sectional  area  of  the  fiber.  The  result  is 


£  elk":  2 ikm  a‘„(z)  H>m  Ex„dA  +  an* o'  k2  $  H>m(UlnVdA 

n 

+  «o  '  $  m,  £j-(U E„)*  +  dA 

+  I  *„£o  1  §  H'm  e‘V'(UE ny-dA  =  0 


£  A  ■  2 ikma’Jz)  <fi  E>n  HxdA  +  am€ 0“'  k2  $  <UE J'H'dA 

m 


+  am*0 ' 


§> 


e'k”:(UEm):HxdA  -  0. 


These  equations  can  be  recast  to  bring  out  the  concept  of  mode  coupling. 

Let  m  =  1  in  Eqs.  (21a)  and  (21b).  Then,  adding  the  two  together,  one  obtains 
e*':2  ikia;(z)/,T  =  -IefVaJ*o',*2C+«o,/.?l 

n 

<f  H\  ( U  E  „ )  xdA  +  <f  (UE,  )'HxdA, 


ll2)  =  £  H\  ii(UE,)'+r^-(UE „V  dA  +£  -4-(UE,)'  +  ^-(UE,)1  HxdA.  (22c) 
]n  *  0x2  9x9.y  y  9  y2 


Define 


Am(z)  =  e,km'am(z)  and  AL(z)  =  a'm(z)e +  ikm  A„,(: 
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Since 


a{{z  )c 


tk„,: 


Z  /•‘fl,fco,*2/|'i’  +  «0 

2  ikiHf 


Cl 


(24) 


=  i  Z  *u  A« 

n 

it  is  seen  that 


Al ’<->  =  ik\  Al(z)  +  '  Z  *l»  An 

n 


(25) 


and 


«0 


'k2l" 


u  +  *0 


1  ;(2) 
*  1« 


2k 


I  Ml 


(26) 


Similar  formulas  are  obtained  for  m  =  2.  The  quantities  Km„  are  the  mode  coupling  coefficients.  They 
describe  how  the  dielectric  polarization  caused  by  the  dyadic  U  makes  the  expansion  coefficient  am  of 
one  mode  depend  on  the  coefficients  of  all  other  modes.  A  particular  case  of  mode  coupling  of  primary 
interest  here  is  the  photoelastic  effect,  in  which  U  describes  the  interaction  of  light  and  mechanical 
stress.  This  is  described  next  for  the  case  of  torsional  stress. 

Dielectric  Polarization  by  Torsional  Stress 


The  dyadic  operation  UE  yields  a  vector  with  three  components: 


(UF.)'=  Uxx  Ex  +  Uxy  Ey  +  Ux:  E:, 

(27a) 

(UE)  ‘  =  UyxE'+  Uyy  Ey  +  Uyz  E:, 

(27b) 

( U  E ) :  =  U.x  Ex  +  U:y  Ey  +  U::  E:. 

(27c) 

When  the  fiber  is  strained  in  pure  torsion  (Fig.  5),  the  displacements  for  a  twist  r  (angle  per  unit 
length)  are 

ux  =  -rzy  and  uv  =  tzx.  (28) 


The  strains  are  then 


Sx:  =  SiX 


dux  du: 
dz  dx 


=  -ry 


(29a) 


and 


SK  =  ^ 


3  uy  3  u: 
dz  +  3.y 


(29b) 


Now  the  components  of  the  dyadic  U  are  proportional  to  the  strain  components, 


U„  =  gS,r 


(30) 
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and 


/  Vm€u 


/  v  ,  ,  v  ~V/M*o  ,  ,  k 

-*t  V”- - y0(r) - 7 —  yo(r) 

A,  r  M 


“*T  *7 


7„(r)y„(r) 


</.4. 


(32d) 


We  omit  for  the  moment  the  calculation  of  the  third  terms  in  (21a)  and  (21b).  Adding  Eqs.  (32a)  and 
(32c),  one  obtains 


/|J'  =  -  ■*?  HT  2ir  f  r2Jt)(r)jit(r)dr. 

k  ,  JJ  0 


A  |  /j. 

Integration  by  parts  leads  to 

. R 


(33a) 


X  Uti(r)  j{t(r)]r2  dr  = 


Jn2(r) 


l 


K  y02  (r) 


2rdr 


Hence, 


“  “X  Jn(r)rdr. 


'^~r  gr  2-rr  f  Jq  (r  )rdr. 


A|M 


Next  we  calculate 

$  H\  E\  dA  =  In  J  R  J$  (r)rdr, 

§>  E\  H]  dA  =  0. 

<f  H\  E)  dA  =  0, 


(33b) 

(34a) 

(34b) 

(34c) 


and 

J  E'2  H)  dA  =  -  In  f*  J 02  ( r )  rdr. 

We  combine  (34a)  and  (34b)  and  combine  (34c)  and  (34d),  and  get 
I  IT  ”  2n  f  ——7-^  Jq  (r)rdr 

•'0  11 

and 


'22  “  “  1 11 

These  quantities  will  now  be  used  to  find  a'm  of  am. 

As  an  example,  let  n  -  2.  Then 

«n  1  «-?  r  gr2n  f  J$  (r)rdr  +  /« »'  _£!_  27r  C  j  (r)j  (r)dr 

K  \  fJi  **  w  ij  k  1  0 


*12  : 


c  I 


4ttA  ,  f  fJ02  (r)rdr 

fX 


(34d) 


(34e) 

(34f) 
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Now  when  r  is  very  small,  Jo(r)  —  0,  hence  we  neglect  /if.  Thus, 


K'2  2  k  \  *T' 


g  =  P4Art$  €„  and  /?„  1  =  A/A,, 


/  |  r  4  *r  , 

*12  ~  T  €«  ~  P44*0  *<>  *  T  ^44  ()  • 

2  /lo  2 

This  result  is  identical  (except  for  sign)  with  that  of  Ulrich  and  Simon  [9] . 

Since  m  and  n  are  dummy  indices  when  they  appear  in  Eqs.  (21a)  and  (21b),  it  is  seen  from 
(340  that 

K  12  **  “*2I* 


We  return  briefly  to  the  third  terms  in  (21a)  and  (21b).  They  both  cancel  for  m  =  1,  n  =  2: 


n  Jk2: 


e‘k 2~  (UE 2y  dA  =  <£  J()(r)  ik2  -^-l—gry  ■  0  +  gTx  y0(r)l 

*  M  d* 

—  —  grik2  - ~ ^  */q  (')  ^  e 'k:' 


'  ^  +  grx  *  0]-^  (~)J0(r)dA  -  /A,  e'*1'  £7  7,?  (r) 


Since  A,  =  A 2  (modes  are  degenerate),  the  z-components  cancel. 

We  next  consider  the  mode  coupling  for  the  case  of  m  =  1,  n  =  1: 

//j  (U  Ej )  v  ^  =  <£)  70  (r)  — ~L  -  grr  sin  <£  —  cos  <£70  (r)  . 

J»  2rr 

sin  <£  cos  <t>  d<t>  —  0,  hence  this  integral  is  zero  and 

0 

(UE,)1  //f  <£4  =§grr  cos2<*>  —■  j0(r)  ■  0=  0, 

# «  &  «*■>'  * -$'*<'>  b  1 1  H  * 

f  !&;■  (UE'>‘  «*  "  -  #  !-*'>■  7  JT  '-''’I"  0  - 0 


^  ~~  gr  -j- J0(r)  dA-  o. 
fx  oxoy  r  k\ 


-gry  ^  70(r)  dA  •  0  =  °- 

o.vflx  r  A  i 


(f)  H\  -~r~  /'"(UE  )Y  dA  =  £  — ^~r~  Jpir )  ik\e'k]‘  [-  gry  J{\(r)  -f  grx  •  0]  dA  —  0 

rivri?  *  n  fix 
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and 

§  ~T-  /'■’  (UE,)--  H\  =  0. 

J  OYOZ 

From  this  it  is  seen  that  km  —  0.  Also,  since  m  and  n  are  dummy  indices  when  used  together,  it 
is  seen  that  k22  0. 

Discussion :  When  m  =  n,  k mn  expresses  the  condition  that  the  energy  density  E  D  is  no  longer  a 
sum  of  diagonal  terms,  because  each  component  of  a  single  mode  is  coupled  to  the  other  two  com¬ 
ponents  through  U  .  The  mode  is  said  to  be  detuned .  When  m  ^  n  the  disturbance  couples  a  com¬ 
ponent  of  one  mode  to  the  components  of  the  second  mode.  Such  coupling  is  forbidden  in  the  undis¬ 
turbed  material,  but  it  is  the  origin  of  experimentally  observed  twist  effects. 


Evolution  of  Electric  Vector  Polarization  in  the  Case  of  Pure  Torsion  of  a  Fiber 


The  set  of  coupled  equations  in  Am  is  now  applied  to  two  degenerate  modes: 

dAx 

-  =  IK  m  A  i  +  /K  j  2  A  ? 

dz 


and 


dA2 

dz 


*  1*2 1 


A  i  +  iK 


A2. 


Since  ,4]  and  A2  are  complex  numbers,  a  convenient  transformation  of  this  can  be  made  in  terms  of 
ii/  and  </>,  representing  their  normalized  amplitude  difference  and  th^ii  m  rmaiized  phase  difference, 
respectively.  The  mean  amplitude  C  and  mean  phase  do  not  contribute  to  mode  coupling  and  so  are 
not  considered.  By  further  introducing  the  spatial  beat  frequency  25  =  -  Ki2,  and  assuming  a  loss¬ 
less  fiber,  k  12  =  kJi  =  Ulrich  [10]  shows  that  the  angular  velocity  11  and  direction  of  the  axis 

R  of  rotation  are  given  by 

n  =  2(82  +  k2>,/2  =  |(k„  -  k22)2  +  4k)jk  I |/2. 

„  8  *11  _  K22 
24/ R  =  arctan  —  —  a  re  tan  — =====  , 

K  74*12*21 

and 

2<t>  R  —  arg  12  =  y . 


Thus  an  evaluation  of  the  mode  coupling  coefficients  k mn  gives  directly  the  location  of  the  axis  of  rota¬ 
tion  and  the  rotation  per  unit  of  length. 

These  results  may  now  be  used  to  find  the  magnitude  and  direction  of  the  axis  of  rotation  in  the 
case  of  pure  torsion: 

In  |  -  |4  k,2k  J, 1 1/2  588  2  kI2|  =  r  \p44\n$  (units:  rad/m), 

2 x  =  2&  r  =  arctan(O)  =  0, 


and 

24/  *  2<f>R  =  arg  *  l2  =  -90° 

Since  \<t>  R\  =  45°,  it  is  seen  from  the  properties  of  the  Poincare  sphere  S  that  the  two  modes  are  in 
time-phase  quadrature ,  hence  the  polarization  is  circularly  birefringent.  Also,  since  4t  R  =  0,  it  is  seen 
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Periodic  Dielectric  Polarization 

When  the  dielectric  polarization  is  periodic  in  space  and  time  it  can  be  represented  in  real  form  by 


and 


P=  UE=  U(r,  /)EcosO 


1  V 


(35a) 


(35b) 


in  which  K  is  the  spatial  wavenumber  of  the  elements  of  U  and  ft  is  the  associated  angular  frequency. 
With  this  form  of  P  the  propagation  of  E  is  found  from  Eq.  (7)  to  be 


V2E  -n't 


82E 

0  8, 2 


82 

U  (r,  /)E  cos  D 

/  -  \ 

dt2 

V 

(36) 


+  ««'V 


U  (r,  /)E  cos  D 

z 

1 - 77 

V 

In  the  case  of  a  homogeneous  glass  optical  fiber  under  mechanical  stress  an  approximate  model  for  E 
can  be  obtained  if  one  assumes  that 


and 


Then 


V«0  —  0, 

E  -  E(r)  e ,  m  »  ft, 
m'  ~  Mo. 

U  (r,  /)  — »  U  (r). 


V2E  (r)  +  ArE  (r)  =  — «0  V'lU  (r)E]  ft2  cos  ft 


l  - 


(37a) 

(37b) 

(37c) 

(37d) 

(38) 


+  «o'  V 


U  (r)E  cos  H 

t  -  \ 

V 

From  this  it  is  seen  that  all  electric  field  polarization  effects  discussed  earlier  in  this  report  will  become 
periodic  with  period  T  =  2 rr/Sl. 

ANALYSIS  BY  PHOTOELASTICITY  OF  ELECTRIC  FIELD  POLARIZATION  OF 
SINGLE  MODE  FIBERS  UNDER  MECHANICAL  STRESS 

In  previous  sections  of  this  report  the  optical  fiber  was  assumed  to  be  carrying  two  degenerate 
modes  although  its  physical  construction  was  "single  mode."  Also,  the  dielectric  polarization  UE  was 
assumed  to  be  caused  by  applied  elastic  strains.  We  next  consider  the  case  of  single-mode  fibers  carry¬ 
ing  only  one  mode  and  assume  UE  to  be  caused  by  applied  elastic  stresses.  The  theory  of  the  electric 
field  polarization  in  this  case  is  the  classical  one  of  photoelaslicity.  Our  objective  in  the  following  sec¬ 
tions  is  to  make  numerical  estimates  of  the  changes  in  index  of  refraction  of  optical  fibers  with  applied 
mechanical  stress.  To  have  a  substantial  basis  for  assuring  the  validity  of  these  estimates,  it  is  neces¬ 
sary  to  review  the  fundamentals  of  the  theory  of  photoelasticity. 
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Stress  Analysis 

Mechanical  stress  in  an  elastic  body  is  represented  by  a  matrix  of  tbrces  acting  on  an  elementary 
volume  centered  at  a  point.  The  matrix  T  in  three-dimensions  (a.  y,  z)  has  nine  components  Ttf  (/  ~ 
v,  v.  z\  j  =  a,  y,  z).  The  first  subscript  indicates  the  x,  y,  or  z  component  of  the  force  through  the 
point,  the  second  subscript  indicates  the  coordinate  plane  through  the  point  on  which  the  component 
acts.  In  two  dimensions  (say  coordinates  a;  y)  the  matrix  has  four  components  oi  plane  stress.  In  iso¬ 
tropic  materials  the  off-diagonal  components  are  equal,  thus  reducing  the  state  of  plane  stress  to  only 
three  distinct  components,  namely,  two  normal  stresses,  crA  and  crk.  and  one  shear  stress  rxy(~T\x). 
Through  the  center  of  the  stress  volume,  in  the  plane  of  plane  stress,  there  is  always  one  direction  (the 
trace  of  a  normal  plane)  in  which  the  normal  stress  (/>)  is  a  maximum,  and  an  associated  orthogonal 
direction  (or  plane)  in  which  a  second  normal  stress  (q)  is  a  minimum  (Fig.  6).  These  are  the  principal 
directions  (or  planes)  of  stress  at  the  point.  In  these  directions  the  shear  stress  vanishes.  Rotation  of 
these  directions  about  the  center  causes  the  shear  to  grow  until  it  reaches  a  maximum  at  a  rotation  of 
45°.  When  the  stress  at  a  point  is  three  dimensional,  one  can  create  a  condition  of  plane  stress  by 
choosing  any  direction  through  the  point  and  calculating  the  stresses  on  a  plane  perpendicular  to  this 
direction.  On  this  plane  one  can  then  locate  the  orthogonal  directions  of  the  secondary  principal  stresses 
by  formulas  given  in  later  sections. 

The  existence  of  primary  principal  stresses  in  two-dimensional  stress  and  secondary  principal 
stresses  in  three-dimensional  stress  determines  the  state  of  polarisation  of  light  propagating  through  the 
point  in  a  stressed  body.  This  is  considered  next. 


Propagation  of  Polarized  Light  in  Doubly  Refracting  Materials 

Every  elementary  volume  centered  at  a  point  of  an  elastic  body  under  mechanical  stress  exhibits 
temporary  double  refraction  to  polarized  light,  provided  certain  conditions  are  fulfilled.  Upon  entering 
the  volume  in  a  given  direction  through  its  center,  the  light  ray  may  encounter  a  state  of  local  plane 
stress  on  a  plane  through  the  center  transverse  to  the  incident  direction.  These  are  the  stresses  that 
alter  the  state  of  polarization  of  the  light.  No  other  stress  components  have  any  effect  on  polarization.  The 
light  vector  oscillating  transverse  to  the  direction  of  propagation  sees  only  two  orthogonal  planes  upon 
which  it  can  vibrate,  all  other  planes  being  effectively  prohibited.  These  are  the  planes  of  primary  or 
secondary  principal  stresses,  whose  traces  are  in  the  transverse  plane.  The  light  therefore  splits  into 
two  components:  one  component  travels  with  its  light  vector  parallel  to  one  principal  plane,  and  the 
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second  component  travels  with  its  light  vector  parallel  to  the  orthogonal  principal  plane.  The  speed  of 
travel  of  each  component  is  determined  by  the  value  of  the  principal  stress  it  encounters  along  its  path. 
Since  the  principal  stresses  are  generally  different  in  sign  and  magnitude,  one  component  travels  on  a 
fast  axis  (or  plane)  and  the  second  component  on  a  slow  axis  (or  plane).  If  the  incident  light  is  of  a 
single  frequency,  this  difference  in  speed  between  the  two  beams  displays  itself  as  a  difference  in  phase 
( A | JS(  -  kskm)z  at  each  coordinate  distance  z  of  the  path,  k  being  the  wavenumber  of  the  propagating 
light.  The  relation  of  altered  wavenumbers  A)as(  and  A\lo*  to  the  mechanical  stresses  p  and  q  forms  the 
basis  of  the  theory  of  photoelasticity.  This  is  discussed  next. 

Repeated  experiments  have  established  the  stress-optic  law  that  A(asl  -  Asj(m  is  proportional  to 
p  -  q\  that  is,  the  difference  in  wavenumber  is  directly  proportional  to  the  difference  of  principal 
stresses  in  the  transverse  plane.  If  the  incident  polarized  light  of  single  wavelength  \0  propagates 
through  a  thickness  d  of  a  stressed  elastic  body,  the  two  components  exit  with  a  phase  difference  of 

-  k{)C(p  -  q)d, 

where  A0—  27r/A0,  C  is  the  stress-optic  constant  of  the  material,  and  p  —  q  is  the  difference  of  primary 
or  secondary  principal  stresses  in  the  transverse  plane.  An  estimation  of  this  phase  difference  and  the 
state  of  exit  polarization  of  the  two  component  beams  requires  careful  consideration.  This  is  now 
reviewed. 

Polarization  Effects  in  Stressed  Elastic  Bodies 

Assume  first  that  the  incident  light  is  plane  polarized  at  an  angle  a  with  the  x  axis  in  an  xy  coordi¬ 
nate  system,  and  assume  that  the  coordinate  axes  are  principal  axes  of  stress.  The  plane  xy  is 
transverse  to  the  incident  ray  (Fig.  7).  The  incident  light  vector  is  represented  as  a  simple  sinusoid, 

s  =  a  cos  <»t. 


It  splits  into  two  unique  components  (no  other  components  effectively  propagating)  initially  synchron¬ 
ized  (that  is,  having  the  same  temporal  phase), 

sx  =  a  cos  a  cos  wt  and  sv  —  a  sin  a  cos  tot. 

These  components  travel  along  the  principal  planes  of  (two-dimensional)  stress  at  different  speeds,  Vx 
and  Vv.  After  a  propagation  distance  d  they  become  unsynchronized. 
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a\  =  a*  cos  (t ot  +  A</>  | ) , 
—  b*  COS  {cot  + 

a*  =  a  cos  a  , 
b*  =  a  sin  a , 

OJ 


A*,  - 

A</>2  = 


K," 

CO 


1 


( A<^> ,  -  A <t>2)/2  =  A<£{/  =  k o  C(p  —  q)d. 


In  the  theory  of  photoelasticity  the  angle  A <t>d  is  measured  to  find  the  stress  difference  p  -  q.  This  is 
done  by  causing  sx  and  sv  to  interfere  in  a  selected  plane  of  a  polariscope.  Choosing  plane  EF  perpen¬ 
dicular  to  the  incident  light  vector,  one  finds  the  amplitude  of  the  resultant  vibration  to  be  Ill,  p.  124] 

B  =  a  sin  2a  sin  A$j. 


If  p  -  q  is  zero,  B  is  zero,  meaning  that  the  light  transmitted  through  point  O  is  extinguished  when 
viewed  on  plane  EF  in  a  polariscope.  A  similar  extinguishment  occurs  when 

&<t>j=k0C(p-q)d=n7T,  n  =  1,  2.... 

If  the  stress  difference  p  —  q  is  a  function  of  time, 

p  -  q  =  pit)  -  qit) , 


then 


B  =  (a  sin  2a)sin{/c0  C[p  it)  -  qit)]d)f 


i.e.,  the  light  viewed  in  the  polariscope  on  plane  EF  varies  in  amplitude  with  time. 

If  the  propagation  direction  is  z,  then  the  light  vector  of  Fig.  7  at  any  distance  z  =  const,  is  the 
vector  sum  of  two  unsynchronized  components  of  different  amplitude.  Thus,  in  the  plane  z  =  const 
the  total  light  vector  rotates  in  time  about  the  center  in  an  elliptical  motion.  The  stress  difference 
p  -  q  can  then  be  said  to  unsynchronize  the  initially  synchronized  components  of  the  plane-polarized 
light  along  the  principal  axes.  Since  the  unsynchronization  varies  with  distance  of  propagation  dy  the 
major  and  minor  axes  of  the  elliptical  motion  of  the  light  vector  change  continuously  with  increase  in 
path.  However,  at  any  station  z  =  const  a  polariscope  will  find  a  plane  EF  along  which  the  transmission 
will  be  extinguished  whenever  p  -  q  reaches  a  level  such  that  A <t>(i  -  nn ,  n  =  1 ,  2  ...  . 

Next  we  consider  the  case  where  the  incident  light  ray  at  z  =  0  has  two  equal-amplitude  orthogo¬ 
nal  components  which  are  temporally  unsynchronized  by  90°: 

sx  =  a  cos  to  t 

and 


sv  =  a  sin  a)t. 

The  light  vector  rotates  in  a  circle  about  the  center  point  of  incidence  in  the  plane  z  =  const  Assume 
that  it  rotates  counterclockwise.  Upon  propagating  a  distance  d  through  a  stressed  elastic  body  in  which 
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P  q  m  the  transverse  plane  is  finite,  the  two  components  of  the  light  vector  are  further  unsynchron- 
l/ed  by  the  temporal  phase  A<£(/  defined  above.  Assume  a  reference  axis  y'  along  which  the  incremen¬ 
tal  phase  vanishes.  Then 

s'  =  a  cos(tu/  4-  =  a  sin (co /  +  it/ 2  + 

and 

si  =  a  sin  cot. 

It  can  be  shown  [12,  p.  94]  that  the  motion  of  the  resultant  light  vector  at  a  local  station  where 
r  —  const  can  be  resolved  into  two  oppositely  rotating  circular  motions: 

si ' 1  =  A  cos  (cof  +  y),  s/1*  =  A  sin  (wf  +  y) 

and 

\<2)  =  —B  cos  (t ot  4-  ft),  si2'  =  B  sin  (o»/  4-  ft), 

where 


2 

4*  *  ~y  1  +  sin  y  4-  <t> j  , 


and 

cos  y  +  <t>ti 

tan  ft  —  - - r . 

1  -  sin  y  +  (bj 

The  combination  uf  the  two  motions  is  equivalent  to  elliptical  motion.  Thus,  the  components  of 
incident  circularly  polarized  light  are  temporally  unsynchronized  more  or  less  than  90°  by  mechanical 
stresses  p  -  q,  making  the  resultant  elliptically  polarized.  As  before,  the  angular  location  of  the  major 
and  minor  axes  of  the  traced  ellipse  varies  continuously  with  distance  of  propagation.  Several  possible 
conditions  of  polarization  can  be  obtained: 

•  If  the  phase  difference  corresponds  to  an  integral  number  of  wavelengths  (  2mr .  n  =  1, 
2  . . .),  no  change  in  the  initial  circularly  polarized  light  results. 

•  If  the  phase  difference  corresponds  to  a  half  wavelength  (plus  any  number  of  integral 
wavelengths),  then  A  =  0,  B  =  a,  and  the  initial  circularly  by  polarized  light  emerges  with 
reversed  direction  of  rotation. 

•  If  the  phase  difference  is  arbitrary,  A  and  B  are  finite  and  different,  and  the  initial  circularly 
polarized  light  emerges  elliptically  polarized. 

Finally,  we  consider  the  case  where  the  incident  light  has  two  orthogonal  components  with 
unequal  amplitudes,  temporally  unsynchronized  by  90°  so  that  the  resultant  light  vector  traces  an  ellipse 
in  a  plane  normal  to  the  optic  ray: 
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and 

If  one  sets  [12,  p.  94] 


and 


then  it  will  be  seen  that 


sx  =  a  cos  oji 


sv  =  b  sin  oji. 


il 

*1- 

a2  +  b2  +  lab  sin 

f 

il 

a2  +  b2  —  lab  sin 

f  +<*,/ 

a  cos 


tan  y  —  — 


TT 

J  +<t>u 


b  +  a  sin 


f  +*« 


a  cos 


tan  5  = 


b  —  a  sin  —  -f  <t>  d 


•  If  <t>j  corresponds  to  an  integral  number  of  wavelengths,  no  change  in  polarization  is  to  be 
expected. 

•  If  <t>d  corresponds  to  a  half-wavelength,  then  the  initial  elliptical  polarization  emerges  ellipti- 
cally  polarized  with  reversed  rotation. 

•  If  <t>d  is  arbitrary,  the  emerging  light  is  also  elliptically  polarized  with  different  sizes  and  direc¬ 
tions  of  major  and  minor  axes. 


Summary 

This  completes  our  brief  review  of  aspects  of  photoelasticity  which  are  relevant  to  the  study  of  the 
effects  of  mechanical  stress  on  polarization  of  incident  light.  Photoelasticians  use  their  tools  to  measure 
stresses  p  and  q.  The  same  tools  in  optical-fiber  studies  require  the  stresses  to  be  known.  We  assume 
in  subsequent  sections  that  the  stresses  p  and  q  can  be  estimated  in  specific  examples.  When  the 
stresses  are  evaluated  we  can  make  numerical  estimates  of  phase  changes  A^.  These  are  also  esti¬ 
mates  of  changes  in  index  of  refraction  A n  caused  by  stress  birefringence.  This  is  explained  later  in  the 
report. 

APPLICATIONS 

Case  1.  Applied  Stress  on  the  Fiber  is  Axisymmetric 
and  Everywhere  Uniform  Along  the  Optical  Path 

Here  the  applied  stress  is  a  scalar  pressure  in  the  medium  which  exerts  a  force  normal  to  the  sur¬ 
face.  The  radial  stress  cr ,  and  the  tangential  stress  on  a  polar  elementary  volume  are  equal  to  each 
other  and  equal  to  the  applied  pressure  [11,  p.  55].  The  principal  stress  difference  p  -  q  is  everywhere 
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zero.  Hence  there  is  no  birefringent  effect  for  an  optic  ray  on  the  axis  of  the  fiber.  Thus,  acoustic 
pressures  in  the  sonar  range  of  frequencies  do  not  alter  the  state  of  polarization  of  light  propagation  in 
them.  Similarly,  hydrostatic  pressures  do  not  affect  the  phase  of  the  light  beam.  It  also  does  not 
matter  what  the  time  variation  of  the  applied  stress  is,  since  no  change  in  polarization  of  the  incident 
light  will  occur  as  long  as  the  stress  is  axisymmetric  on  the  fiber. 

Case  2 .  Applied  Stress  on  the  Fiber  is  Axisymmetric  but  Varies  Along  the  Optic  Path 

Here  again,  as  long  as  the  applied  stress  is  axisymmetric  (that  is,  an  applied  scalar  pressure),  there 
is  no  birefringent  effect  for  axial  optic  rays  regardless  of  the  nature  of  the  variation  of  stress  along  the 
path  of  propagation. 

Case  3.  Tensile  or  Compressive  Stresses  are  Applied  to  the  Ends  of  the  Fiber 

In  this  case  the  fiber  is  uniformly  stretched  or  compressed.  All  stresses  inside  the  fiber  on  cross 
sections  normal  to  the  fiber  axis  are  parallel  to  the  propagating  optic  ray.  Since  there  are  no  transverse 
components  of  stress  for  an  axial  optic  ray,  there  is  no  birefringent  effect.  When,  however,  the 
incident  optic  ray  is  oblique  to  a  selected  plane  of  stresses,  there  is  an  effect.  This  is  discussed  next. 

Case  4.  Oblique  Incidence  of  Optic  Ray  and  Pure  Bending 

Figure  8  shows  stress  components  on  a  three-dimensional  elementary  volume  centered  at 
xOt,  y,  z)  relative  to  the  origin.  The  secondary  principal  stresses  at  x  in  direction  i  are  defined  (as 
noted  previously)  as  *he  principal  stresses  resulting  from  stress  components  which  lie  in  a  plane 
transverse  to  i  through  x.  These  are  designated  (p' ,  q As  an  example  of  secondary  principal 
stresses  choose  the  z  direction.  When  a  A  and  cr  v  are  present,  one  finds  [13,  p.  90] 


Pig  K  —  Stress  components  on  a  three-dimensional  elementary  volume 
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The  orientation  of  p '  and  q *  at  center  point  x  on  the  xy  plane  is  along  orthogonal  axes  rotated  through 
an  angle  «  relative  to  the  x  axis. 


(tan  2e ).  - 


2ru 


An  optic  ray  propagating  in  the  coordinate  direction  z  obeys  the  following  rules:  (l)  Stress  components 
subscripted  with  the  letter  z  have  no  effect  on  the  polarization.  Thus,  &  z  ,  t  v:,  and  t  x:  have  no  effect. 
In  particular,  pure  torsional  moments  Mz  which  produce  rv.  and  r  x:  have  no  effect.  12 1  If  the  orienta¬ 
tion  of  the  secondary  principal  stresses  rotates  as  the  optic  ray  advances,  the  components  of  polarization 
also  rotate. 

When  the  optic  ray  is  obliquely  incident  on  the  elementary  volume,  the  birefringent  retardation  is 
calculated  from  transverse  components.  A  simple  example  is  a  bar  stretched  by  pure  tension  in  the  x 
direction  (Fig.  9)  and  pierced  by  an  optic  ray  in  the  z  direction  parallel  to  the  xz  plane.  Since  stress 
components  subscripted  with  the  letter  z'have  no  effect  on  polarization,  the  secondary  principal  stresses 
are 

p'  -  <J t  sin2  £ 
and 

Q* =  0. 


Defining  retardation  Rd  =  m\  as  the  number  m  of  wavelengths  X  (generally  not  an  integer)  of 
difference  in  the  phases  of  the  birefringent  components,  one  finds  that  the  retardation  along  oblique 
path  /'  is  then  [14,  p.  339], 

Rd  =  mk  =  Ct  'cr ,  sin2  £  =  Ct<r ,  sin  £ , 

where 


t  *  t'  sin  £ . 

This  formula  will  be  used  in  later  numerical  calculations.  We  next  apply  oblique  incidence  to  pure 
bending  of  a  bar  lying  along  the  x  direction  and  bent  by  pure  moments  Mx  at  a  cross  section  (Fig.  10). 
Let  the  optic  ray  be  incident  in  the  plane  of  the  stresses  cr  at  an  angle  £  with  the  x  axis.  Assume  that  it 
enters  an  elementary  volume  at  y  and  exits  at  y0.  The  secondary  principal  stresses  at  the  center  of  the 
volume  are 
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q'  -  0, 


in  which  L  is  the  areal  moment  of  inertia  of  the  cross  section  about  axis  z.  The  retardation  is  [14,  p. 
348] 


Rd  =  m\ 


c  L  w-w- 


CM, 


sin  £ 


(y2  -  W ) 
2 


Case  5.  Fiber  Subject  to  Both  Bending  and  Shear 

A  fiber  of  circular  cross  section  is  subject  to  a  bending  moment  Mx  and  a  shear  Vv  at  the  ends  of 
a  segment  whose  thickness  is  in  the  x  direction  (Fig.  11).  The  normal  stress  at  distance  y  from  the 
neutral  axis  is 

Mxy 


The  shear  stresses  t  for  all  points  on  line  pp  at  a  distance  y\  from  the  neutral  axis  all  point  toward  the 
intersection  5  of  tangent  lines  pS.  Each  r  has  two  components  [14,  p.  121], 


rxy  =»  r  cos  P 


V(R2 —  yf) 
31. 


and 


rKZ  =  r  sin  p 


3/. 


where 


VR^R2  - 
31. 


y  i 


Jr 


+  t£. 
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The  principal  stresses  lie  in  a  plane  normal  to  x  =  const.,  rotated  through  an  angle  c i>  about  the  r  axis, 
where 

tan  2 <t>  =  — . 

<*x 


and 

(T 


To  summarize  the  effect  of  these  stresses  on  the  polarization  state  of  the  optic  ray:  (1)  An  optic  ra> 
propagating  exactly  in  the  a  direction  normal  to  the  cross  section  experiences  no  change  in  polarization, 
since  all  stress  components  have  the  letter  v  in  their  subscripts.  (2)  If  the  optic  ray  is  parallel  to  the.tr 
plane  at  an  angle  {  to  the  x  axis,  the  secondary  stresses  are 

p  —  q'  —  (p  -  </)sin:£. 

The  retardation  along  ray  path  /'  is  then 

mk  =  dp'  —  q')t\ 

where  p\  q\  and  /'  are  functions  of  the  distance  y  from  the  neutral  axis.  If  the  optic  ray  points  in  a 
direction  fixed  by  direction  cosines  ar,  and  y,  (relative  to  a,  y,  and  r)  and  the  stresses  p  and  q  lie  in 
a  plane  through  the  stress  point  normal  to  the  at  plane,  the  direction  cosines  of  />  being  a,,.  f3r,  and  0, 
and  those  of  q  being  fiir  and  0,  then  the  components  of  p  and  q  transverse  to  the  optic  ray  are 
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cos  <b  =  cos  a,  cos  a„  +  cos  ft,  cos  fi., 


cos  f  =  COS  tt,  cos  +  cos  ft,  cos 


I  he  retardation  is  then  given  by 


m\  —  C(/>  sin'  </j  —  q  sirr  Of' 


C  ase  b.  fiber  Wrapped  on  a  Mandrel  (151 

A  circular  fiber  of  diameter  d  is  wrapped  on  a  mandrel  of  radius  r.  A  tangential  force  T  is  applied 
to  the  two  open  ends  of  the  fiber  to  keep  it  tight.  The  stresses  on  a  cross  section  at  angle  4>  are 
analyzed  with  relation  to  the  forces  N  and  moments  M  shown  in  Fig.  12.  These  are 

M  =  Tr  ( 1  —  cos  4 > ) , 

A  =  T  cos  </> . 


Q  =  T  sin  4). 

A  cross  section  of  the  ring  exhibits  shear  and  bending  (Fig.  13).  The  neutral  axis  of  bending  is  orthog¬ 
onal  to  the  plane  of  the  ring  and  passes  through  the  axis  of  the  fiber. 


12  —  Forces  on  a  fiber  wrapped 
on  a  mandrel 


Shear  (Jand  bending  moment  M  on  a 
cross  section  &  =  const 
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Let  v  represent  the  coordinate  axis  of  the  fiber,  and  y  the  perpendicular  distance  of  any  point  from 
the  neutral  axis.  Then  the  bending  stress  is 

My 


A  ' 


where  A  -  and  M  =  2  Tr  sin:  .  Thus 
64  2 


128  Try  sin2  ^ 

=  - ^ - • 

The  shear  stress  consists  of  two  components  (see  previous  section), 
cl2 


(A 


—  V 


3  A 


64  J_ 
3  dA 


d2 


A  sin  </> 


and 


Qy 


1/2 

d2 

T  -  -‘"1 

64 

T  ~  -r 

1/2 


3  A 


3  ' 


A  sin  <l> , 


and 


/  T-  ~  32  1 

-  >/r;,  +  ra  -  T  p- 


l  “ 


1/2 


A  sin  <£. 


The  principal  stresses  at  point  (y,  z)  lie  in  a  plane  normal  to  the  .yz  plane: 


and 


Q  - 


2 

+  T2 

2 

c rx  1 

2 

+  T2 

2  , 

1/2 


1/2 


They  are  on  two  orthogonal  axes  rotated  about  the  z  axis  through  an  angle  «  given  by 


tan2«  =  — — . 

cr , 


In  terms  of  component  shears  we  can  write  component  principal  stresses  defined  by 


(P-  2 


64>Tr  sin 


+ 


and 


<P  -  «/>.-,  =  2 


64 y7>  sim 


32 

3 


32 


A  sin 


</2 


-  y i 


1/2 


A  sin 


</2 

T  - 


1/  2 


r/4 


1/2 
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The  effect  of  these  stress  components  on  polarization  depends  on  the  angle  of  incidence  of  ,  e  optic 
ray.  If  the  ray  is  along  x  exactly  normal  to  the  cross  section  the  effect  vanishes,  since  the  subscripts  of 
the  stress  components  contain  the  letter  x.  Numerical  calculations  for  optic  rays  obliquely  incident  on  a 
cross  section  are  given  in  a  later  section. 

Case  7.  Fiber  Subject  to  Random  Stresses 

Let  the  stress  fluctuations  (defined  as  8  s  =  8  (p  -  q))  be  randomly  distributed  along  path  length 
x .  Assume  that  these  stresses  can  be  represented  by  a  single  scalar  function  of  x.  Then  the  mean- 
square  phase  after  length  L  is  reached  is 

<h<f>j>  =  <  [ J*0  kQC(x)bs (x)dbtj  > 

=  kQ  C2<|j*o  8s(x)dx]  > 

if  C(x)  =  C  —  constant.  In  forming  the  square  on  the  right-hand  side,  we  must  allow  cross  products  at 
various  values  of  x,  x\  etc.  To  evaluate  this  integral  we  use  Taylor’s  method  [16].  By  definition  of  the 
normalized  autocorrelation  function  p, 

f  <8s  (x ) 8s  (x')  >  dx' =*  <8s2>  f  p(x  —  x’)dx'. 

J  o  0 

Since  we  always  take  the  correlation  p  to  be  an  even  function,  we  can  write  the  right  hand  side  as 

<Ss2>  J*p(g)d(. 

Now, 

f  <8s  (x)  8s  0c')  >  dx'  —  <8s(x)  f  8  s  (x ')  dx '  =  <8sGc)^>, 

JO  JQ 

where 

X  =  f  8  s  (.v ')  dx  \ 

J  o 

But 

<8s(x)jt>  =  L-r  <*2> 

2  dx 

Thus,  over  a  length  L, 

<X2>  =  2<8s2>  /o  J*p(£)d{d> k, 

hence, 

=  A:02  C2  ■  2<8s2>  p  (£ )  dt  dx. 

Two  important  cases  can  be  distinguished: 

(a)  Suppose  path  length  L  is  short  enough  so  that  p(£)  is  essentially  unity  over  the  range  of 
integration.  This  is  the  coherent  case.  Then 

<&<b  f>  =  koCL  V2 <8s2>  , 

that  is,  the  rms  value  of  the  phase  is  proportional  to  the  length  of  the  fiber. 

(b)  Assume  L  is  so  long  that  the  stress  at  the  end  of  the  fiber  is  uncorrelated  with  the  stress  at 
the  beginning  of  the  fiber.  This  is  the  incoherent  case.  Then 

J<Ml>  =  k0C  -,/2<8s2>LL(., 
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where 

Lt  =  f>(t)d%  .  L,  «  L. 


In  this  case  the  rms  phase  is  proportional  to  the  square  root  of  the  length  of  the  fiber.  The  symbol  L{ 
represents  the  integral  scale  (or  correlation  length)  of  the  random  field. 

In  both  cases  we  need  the  variance  <&s2> .  This  can  be  obtained  from  the  spectral  density  hSs  of 
the  stress  distribution.  By  Rayleigh’s  theorem, 

f  ^  |8s(x)|‘r/.v  =  L  <8s->  =  2ir  f  |6S,(«)Prf«. 

or 

<8s‘>  =  f  \8S,  (a )  P  da. 

L  J 

The  spectral  density  must  be  determined  by  measurement. 


Case  8.  Fiber  Squeezed  by  Diametral  Forces:  Circular  Disk  Model  111,  p.  82] 


A  unit  thickness  of  fiber  is  modeled  as  a  circular  disk  (Fig.  14).  Let  two  equal  and  opposite  forces 
P  (units:  N/m)  act  along  a  diameter.  The  stress  distribution  is  radial. 


(T  t  — 


2 P  cos  0 
tt  r 


(In  =  0,  Trfi  =  0. 


c.  COS0j  COS  9 2  1 

Since  —  =  — - —  “  ~  and  since  rr9  =  0  the  principal  stresses  on  the  circumference  are 

IP 

cr  min  max  (T  r  ^  • 

7 T  U 


For  equilibrium  one  must  superimpose  tensions  2P/tt  D  on  each  of  these  radial  stresses.  Suppose  next 
the  stresses  on  a  horizontal  diametral  section  are  considered  (Fig.  15).  The  normal  stress  on  plane  ac  is 
contributed  by  two  components  of  the  <rf’ s, 
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On  planes  bc\  dc\  ab ,  and  ad  the  shear  stresses  are  all  zero.  The  stress  in  the  x  direction  is 

4 P  cos  9  ,  2 P 

7 r  D 


( r  v  = - -  sin2#  + 


Thus  o-  v  and  cr  v  are  principal  stresses,  and  their  difference  along  a  horizontal  diameter  is 

4Pcos#  .  2n  ,  2 P  t  4Pcos9  2P 

p  -  q  = - sinz#  4-  — —  +  -  cos^# - — 

nr  7 r  D  nr  n  D 


4P  cos  0  (cos  2q  _  sjn2^)  _  AE  cos  ^(i  -  2sin2#). 
nr  nr 


At  the  center,  x  =  0  =  y. 


P  ~  Q  = 


8  P 
n  D 


At  the  boundary,  x  —  -y ,  y  =  0, 


P  ~  <7  -  0. 


Numerical  calculations  of  this  case  are  given  in  a  later  section. 

Case  9.  Pure  Torsion  Hi,  p.  2351 

A  fiber  of  circular  cross  section  (radius  a)  is  subject  to  a  torsional  moment  M,  (Fig.  16).  On  a 
plane  r  =  const  ,  the  shear  stress  is  resolved  into  two  components. 
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Fig.  16  —  Fiber  subject  lo  pure  torsion 


and 


2Mfy 


2  M(x 


In  polar  coordinates  r,  a , 


2 M.  rsin  a 


it  a 


and 


The  combined  shear  is 


and 


2  M,  r  cos  a 


=  (r  +  r  ;,) 


2  \  1/2  _ 


2M,  r 


7 t  a 


2M, 


16  M, 

t Td} 


The  angle  of  twist  per  unit  of  length  (fl)  is  given  by 

M,  32  M, 

^  GIp  G7rdA 

where  G  is  the  modulus  of  rigidity  and  lp  is  the  polar  moment  of  inertia  of  the  area  of  the  cross  section. 


Assume  that  the  optic  ray  is  incident  obliquely  (at  angle  9)  to  the  cross  sectional  area  parallel  to 
the  xz  plane  (Fig.  17).  Resolve  stress  r  x:  into  components  parallel  to  and  transverse  to  the  ray.  The 
only  effective  component  is  transverse: 

p'  =  (Tx  ==  tx:  sin2  9  ;  </'=0. 
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Next  consider  tvz.  It  has  a  pair  tzx  which  lies  in  the  at  plane  (Fig.  18).  The  transverse  component  of 
t :  v  is  tX}  —  t:}  sintf.  Thus,  the  secondary  principal  stresses  at  point  ( x.y )  of  the  circular  area  of  the 
cross  section  are 

P'  “  y/*rj  +  4  t- 

and 

4'  =  -j-  “  y  s/vJ  +  4t(2,  ; 

hence, 

—  </'  =  2sin  ft  Vt  cos20  4  r  2 

M(  I — r - — = 

=  4sin  - -  vyzco$'t)  4-  xl 

7T  a 

4sin  M f  r  , — ^ 5 - r— 

=  - - —  v sin~a  cos"0  4  cos  a. 

7r  a 

For  fixed  this  stress  difference  is  maximized  at  r  =  a: 

4sin  M,  r-i - 5 - t— 

(p  ~  q’) mix  =  - : —  v<r  sirra  4  cos" a. 

tt  a s 


Fig  18  —  Oblique  incidence  wiih  r n.  component 


S  NANISM 


When  x  —  0,  namely  when  points  are  taken  along  the  y  axis, 

4sin  0  Mt 


P  -  q  = 


y  cos  to 


and 


4sin  0  M{  2Mt 

~  q\\ max  =  - T -  COS  0  =  - T  Sin2  0. 


TTU' 


When  y  =  0,  namely  for  points  along  the  x  axis, 

[(p  q  ^.v^rnax 


4sin0A/,  32sin  0  M( 


7 TCI 


TTd^ 


This  case  is  computed  numerically  in  the  next  section. 


NUMERICAL  CALCULATIONS 


Stress-Optic  Constant  C 

The  theory  of  photoelasticity  specifies  the  value  of  the  stress-optic  coefficient  in  terms  of  a  fringe 
value /defined  as  follows: 


Let  m  be  the  fringe  order  and  mk  the  phase  retardation,  in  units  of  wavelength,  between  the  fast 
and  slow  waves  in  a  birefringent  material.  Then  the  stress-optic  law  requires  that 


mk  =  C(p  —  q)d. 

p  -  q 


Photoelasticians  measure  maximum  shear  stress  rmax  = 
fringe  value  Fby  the  relation 


For  convenience  they  define  a  model 


P  ~~  q 


=  mF, 


from  which  it  is  seen  that 


Thus, 


F  = 


ICd 


or  C  = 


k 

2Fd 


mX  =  —  (p  -  <?). 

Here  F  is  a  material  parameter  whose  value  for  a  given  material  depends  on  the  wavelength  of  the 
incident  light  and  on  the  thickness.  Usually  values  of  F  are  tabulated  for  a  standard  thickness  of  1  in. 
(25.4  mm).  It  is  then  called  the  fringe  value  f  the  units  of  which  are  (lbs/in.2)  x  inch/fringe  (6.895 
kPa  x  2.54  cm/fringe).  For  any  other  thickness  (or  propagation  length), 

/ 

d 

Thus  the  fringe  order  is 


F  = 


(units:  psi/fringe  or  6.895  kPa/fringe). 


...  (P  -  Q )d 

2/ 


(units:  fringe) 


Now  for  every  retardation  (between  slow  and  fast  waves)  of  one  wavelength  the  phase  retardation  is  27? 
rad.  Hence,  the  total  phase  retardation  in  a  fiber  of  length  L  is 
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,  .  -  7 Tip  ~  U)L  t  .  , 

4> j  —  mk  x  2v  =*  - — -  (units:  rad.). 

(For  consistent  units  />  and  </  are  in  psi  [6.895  kPa],  A  is  in  in.  (2.54  cm].) 

The  value  of  /  in  the  case  of  glass  varies  with  the  composition  and  wavelength  Frocht  [12,  p. 
350]  gives  f  —  490  to  1500  psi  •  in.  (86  to  263  kPa  ■  m)  for  a  wavelength  of  0.5461  fim.  Jessup  [17] 
cites  a  value  for  glass  in  the  range  /  =  800  psi  •  in.  (140  kPa  •  m)  upward  . 

A  diff  erent  approach  to  calculating  C  is  adopted  by  physicists.  They  begin  with 

A <f>d  =  mk0  =  ~~  C(p  -  q)d  =  k0C(p  -  q)d. 


From  this  they  define  an  effective  difference  in  index  of  refraction  A  ncf)  between  fast  and  slow'  rays 
caused  by  mechanical  stress, 

A wett  =  ~  y)  (units:  none). 


It  is  thus  seen  that  the  units  of  C  are  the  inverse  of  the  units  of  mechanical  stress,  namely  m7  N~]. 
The  quantity  C  is  related  to  the  classical  stress-optical  coefficients  q„  of  a  material.  If  the  material  is 
isotropic,  then  [6,  p.  705] 

A (c/ 1 1  ~  q\i)(p  -  q)  (units:  none). 


Thus, 


C 


2 


(</n  -  4Ti2> 


units: 


Numerical  values  of  qu  and  q  )2  are  usually  derived  from  a  knowledge  of  the  strain-optic  constants  pu 
of  the  material.  For  stiffness  tensors  C/k  and  compliance  constants  s,,  one  has 

6 

Pu  =  X  Q,k  Q 

k-  1 

and 

6 

</„  =  X  Pik  sjk ■ 
k-  i 

Since  fused  silica  is  elastically  isotropic,  one  finds 


and 


P n  —  Q\i  C\\  +  2Q\2  C\2 


P\2  <7  II  <^12  +  <7  1 2^1 1  +  C12). 


Now,  for  fused  silica,  n 


1.46, 


CM  =  7.85  x  1010  N/m2 
C \2~  1.61  x  1010  N/m2 


Ref.  18, 


and 


P ii  ^  0.121 

Pn  -  0.270 


Ref.  19. 
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Solving  by  Cramer's  rule,  we  obtain: 

</u  =  3.98  x  10  13  m2/N  =  2.74  x  10  9  in.2/lb 
and 

qn  =  2.78  x  10~12  m2/N  =  1.92  x  10  8in.2/lb. 

Thus  we  estimate  C  to  be 

c  =  <2-74  x  10-9  -  i  92  x  i(r8) 

=  -  2.55  x  10~8  in.2/lb  -  -  3.70  x  10~12  m2/N. 

We  check  this  value  by  use  of  a  formula  from  Frocht  [12,  p.  156] 

C  =  ^7  ”  5  t61  ?a!I!  5  ~  2.58  X  10-8  in.2/lb  (3.74  x  10-|2m2/N). 

2/  2  x  1060 

The  magnitude  is  close  enough.  Thus  if  we  retain  the  units  of  X  in  centimeters,  take  /to  be  1060,  and 
disregard  the  negative  sign  we  see  that  our  computed  values  of  qu  and  q]2  are  verified.  An  additional 
check  on  the  value  of  C  is  found  in  [9],  where  C  for  fused  silica  is  given  as  —3.36  x  10" 12  m2/N.  This 
agrees  closely  with  the  values  noted  above  from  other  sources. 

We  proceed  now  with  the  analysis  of  several  cases  of  elastically  stressed  bodies. 

Numerical  Calculations  of  Fiber  on  a  Mandrel  (Case  6) 

A  fiber  80 /im  in  diameter  with  a  cladding  55  nm  in  diameter  and  a  core  4.5  /xm  in  diameter  is 

wound  on  a  mandrel  which  is  1  cm  in  diameter.  From  the  tension  in  the  fiber  T  we  wish  to  determine 

the  difference  in  principal  stresses  in  the  core  and  the  birefringent  effect  per  unit  complete  turn. 

We  begin  with  an  analysis  of  the  loading.  The  relative  areas  are 

A  tot  =  ~  =  j  x  (80  x  10-6)2  =  5.03  x  103  Aim2(7.79  x  lO^6  in.2), 

A  cladding  =  ~[<55  x  lO'6)2  -  (4.5  x  lO’6)2]  =  2.36  X  103Mm2(3.66  x  10“6  in.2), 

A core  -  y(4.5  X  10~6)2  =  15,9  jum2(2.47  x  10~8  in.2), 
and 

^substrate  =  5.03  X  103  -  (2.36  x  103  +  15.9)  =  2.65  x  103  Mm2(4.105  x  10~6  in.2). 

Assume  that  a  tension  of  170  MPa  (25,000  psi)  is  permitted.  The  tensile  force  on  each  area  is 
7"tot  —  25,000  x  7.79  x  10'6=  0.195  lb  -  88  4  g  (0.867  N), 
i  x  i  rr6 

^cadding  “  Jj/  *  0.1948  =  9.15  x  10“2  lb  =  41 .5  g  (0.407  N) , 

T  47  y  10“8 

rcore»  x  X  0.1948  =  6.17  X  10“4  lb  =  0.28  g  (2.74  mN), 
and 

4  ios  y  in-6 

7\ubs.ra,e=  ~ ~~ *  0.1948  =  1.03  x  10'1  lb  =  46.6  g  (0.458  N). 
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We  now  undertake  a  solution  to  the  problem. 

We  know  from  Case  6  that  the  principal  stresses  are  in  the  XV  plane  and  the  XZ  plane: 


(p  ~  <?),,,  =  2 


64y7>  sin2 

2 


and 


( P  ~  =  2 


v  d4 


64yTr  sin2 


nd4 


2 

1/2 

_i_ 

32  T  sin  <f> 

d 2  3) 

2 

t 

3  d4 

l4  II 

2 

+ 

32  T  sin  <t> 

3  y  d4 

d 2  2 

T~y 

1/2 

1/2 


We  consider  ( p  -  q)xy  first.  Two  approaches  appear  feasible  in  estimating  the  polarization  effect  of 
this  quantity.  In  the  first  we  consider  the  fiber  to  be  homogeneous  and  subject  it  to  a  maximum  ten¬ 
sion  7Yot.  We  then  calculate  the  principal  stresses  at  y  =  4.5/2  =  2.25  /urn,  for  d  =  80  /am.  In  the 
second  approach  we  consider  only  the  core  and  subject  it  to  a  tension  Tcorc.  The  calculations  are  as  fol¬ 
lows.  For  component  ( p  -  q)xv:  T  =  Ttot,  d=  80  /am,  y  =  2.25  /am,  and  <t>  =  90°.  Then  d  =  80 
x  I0“h  x  39.37  =  3.15  x  10~3  in.,  y  -  2.25  x  10“6  x  39.37  =  8.86  x  10~5  in.,  and 

U 


(P  -  ?)*>■  “2  T 


64  x  8.86  x  1Q~5  x  1.97  x  10~‘  x  0.5 


7t  (3.15  x  10'3)4 


+ 


32 


1 


(3.15  x  IQ"3)4 


3.15  x  10‘3 


-  (8.85  x  10“5)2 


1/2 


=  2T(3.26x  10u  +  7.22x  I0'°)l/2 


=  3.65  x  106  T psi,  Tinlb(55.5  T GPa,  Tin  kg,  or  5.66  T GPa,  TinN). 

The  secondary  principal  stresses  at  4°  inclination  of  the  optic  ray  are 
ip'  ~  q')x,y  -  3.65  x  106  T  x  sin2(4°) 

=  1.776  x  104  T  psi,  T  in  lb  (270  T  MPa,  T  in  kg,  or  27.5  T  MPa,  Tin  N). 
The  total  loading  is  Ttoj  =  0.1948  lb  (0.0884  kg  or  0.867  N).  Thus, 

(p'  ~  q')x.y  =  1-776  x  104  x  0.1948  =  3.46  x  103  psi  (23.9  MPa). 

Thus  the  difference  An  between  slow  and  fast  components  is 

An  =  C[(p'  -  </').vl]max  =  -2.55  x  10~8  x  3.46  x  103 

=  -8.83  x  10“5. 

Now  assume  <t>  =  180°;  then 

,  ,  ,[  64  x  8.86  x  10'5  xl.97  x  lo1  xl  1  ^ 

*'  ■  ‘  2| - w  (3.15  x  10"  V — r 

and 

=  7.22  x  106  r  psi.  Tin  IbdlOTGPa,  T in  kg,  or  11.2  T GPa,  HnN). 
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The  secondary  principal  stresses  at  4°  incidence  of  the  optic  ray  to  the  fiber  axis  are 

ip'-q')x.y-7. 22  x  106  T  x  sin2(4  °) 

=  3.51  X  104  r  psi.  T  in  lb(  535  T  MPa,  Tin  kg,  or  54  6  7  MPa,  /  in  N> 

Again  let  T  =  7tot  =  0.1948  lb  (0.0884  kg  or  0  867  N)  Then 

( p '  -  </')v>  =  3.51  x  104  x  0  1948  =  6.85  x  U)lpsi  (47  2  MPa) 

Thus,  t  , 

^n  -  -2.55  x  10  8  x  6.85  x  10'  =  -1.75  x  10  4 

Next,  component  (p  —  q)x:  is  seen  to  be  about  the  same  as  ( p  -  q)K Hence,  for  =  180°, 
the  total  A n  is  about  V2  larger  than  for  either  component;  i.e.,  A n  =  -1.75  x  10  4  x  V2  =  -2.5 
x  10-4  at  a  tensile  force  T  ~  ^TOT  which  is  close  to  the  limit  (that  is,  dose  to  the  ultimate  stress  of  the 
glass  fiber,  taken  here  to  be  170  MPa,  or  25,000  psi). 

Next,  we  make  a  calculation  as  if  the  core  alone  is  wrapped  on  the  mandrel. 

Bare  Core  on  a  Mandrel  (Case  6) 

Let  us  consider  only  (p  —  q)xv,  where 

„  _  ,  ,  =  6£  JL 

HP  Q'x.y  imax  ^  ’ 

Using  conventional  units  ( T in  lb,  r  and  d  in  in.), 
r  =  0.5cm  =  0.5/2.54  -  1.968  x  10”'  in. 


and 

d  =  4.5/x m  =  4.5  x  10“6x  39.37=  1.77  x  10  4in., 

thus 

[(  _  .  i  =  64  T  x  1.968  x  1Q~ 1 

UP  9V>lmax  (K77xl0-«)J 

-  7.23  X  10"  T psi,  Tin  lb(  11.0  TPPa,  T in  kg,  or  1.12  TPPa,  TinN). 

The  secondary  stresses  at  4°  inclination  of  the  optic  ray  are 
((/>’- 7.23  x  10"  Tx  sin2  (41 

-  3.51  x  109  T  psi,  Tin  lb  (53.3  TTPa,  7'inkg,or5.44  TTPa,  TinN). 

Thus  the  increment  in  index  of  refraction  is 

An  -  C[(p’  -  q')x,y]m„  -  -  2.55  x  10"8  x3.51  x  109  T 

-  -89.7  T/lb,  7~  in  lb  ( —  1 98  T/  kg,  Tin  kg,  or  20.2  T/  N,  TinN). 

Let  r-  rcore-6.17x  10‘4  lb  (0.280  g  or  28.6 mN).  Then 

An - 89.7  x  6.17  x  10“4=  -5.53  x  10“2. 

Parameter  Variation  for  Calculations  of  An  of  a  Bare  Core  on  a  Mandrel  (Case  6) 

(a)  Let  the  parameter  be  the  angle  £  between  stress  system  and  optic  ray.  The  functional  depen¬ 
dence  of  An  on  £  is  sin2£.  Let  T  =  Tcorc  —  6.17  x  10-4  lb  (0.280  g  or  2.74  mN).  Then  the  relation¬ 
ship  between  £  and  An  is  given  in  Table  1. 
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Table  1 


i 

A  n 

0° 

0 

1° 

-3.44  x  KT5 

2° 

-1.37  x  10'2 

3° 

—  3.095  x  10'2 

4° 

-5.5  x  10'2 

(b)  Let  the  parameter  be  the  tensile  force  T  =  Tcorc  at  the  ends  of  the  fiber.  Assume  f  =  4°  is 
the  angle  of  incidence  of  the  optic  ray.  Then  the  relationship  between  T and  A n  is  given  in  Table  2. 


Table  2 


7  roi 
(lb) 

^ror 

<g> 

^TOl 

(N) 

T 

•  core 

(lb) 

F«ire 

(g) 

T 

1  core 

(N) 

Art 

1  947  x  10  2 

8.83 

8.66  x  10  2 

6.17  x  10~5 

2.80  x  10  2 

2.74  x  10"4 

1 

la 

La 

X 

10“3 

1  947  x  10  - 

8.83  x  10' 1 

8.66  x  10  1 

6.17  x  10'  6 

2.80  x  10  2 

2.74  x  10'5 

-5.5  x 

10  4 

1.947  x  10  4 

8  83  x  10  2 

8.66  x  10-4 

6.17  x  10'7 

2.80  x  10' 4 

2.74  x  10~6 

-5.5  x 

10" 5 

1  947  x  10  " 

8.83  x  10  ' 

8.66  x  10  5 

6  17  x  10  s 

2.80  x  10' 5 

2.74  x  10'7 

-5.5  x 

10  6 

Fiber  Squeezed  by  Diametral  Press  (Case  8) 

The  effective  birefringent  increment  A n  is  obtained  from  Case  8.  We  have 

A n  =  C {p  -  q)  =  —2.55  x  10~8  in.2/lb  x  — 

tt  D  in.2 

=  -6.50  x  10“8 


Suppose  first  that  the  force  Pacts  on  the  core ,  D  —  4.5 /urn  =  1.77  x  10  4  in.  Then 

-6.50  x  10“‘ 


A  n  = 


1.77  x  10" 


P  =  -3.67  x  10'4P. 


To  estimate  P,  assume  that  <rv  is  the  maximum  allowable  stress,  say  170  MPa  (25,000  psi).  Then  at 
the  center  of  the  fiber  disk 

|<r, I  -  170  MPa  -  25,000 

in.  D 

Thus, 

„  n  x  1.77  x  10-4  in.  x  25,  000 lb/in.2  1L/.  /J/wr  \ 

P  -  - - - - - *  2.32  lb/in.  (406  N/m). 

6 

The  differential  A n  due  to  birefringence  is  then 

An  =  -3.67  x  1(T4  x  2.32  =  -8.5  x  1(T4. 


Average  A n  Along  a  Horizontal  Diameter  (Case  8) 

The  stress  p  -  q  can  be  averaged  along  a  horizontal  diameter  by  integration  over  the  polar  angle 
9: 
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/  v  4 p  fnl 4  {cos  0  —  2cos  9  sin3  W )  c/W 

J0  - 7/4 - 

_  2.4  P  =  4.8  P 

7T  r  it  D 

Hence,  the  differential  An  due  to  birefringence  is 

(Aw)av  =  -2.55  x  1(T8  x  — 4  8  x  2.32  =  _5  ,  x  |0  4 

A  it  x  1 .77  x  10  4 


Pure  Torsion  (Case  9) 

Let  us  use  the  last  formula  of  the  analysis  of  Case  9  to  find  An.  Now 


(  )  max 


7T  d* 
16 


Let  d  =  80 Atm  =  3.15  x  10'3  in.,  and  assume 

Tma*=  2.8  X  10-3C  -  2.8  X  10'3  x  4.5  X  106 
=  1.3  x  104 psi  (89.6  MPa). 


Here  we  have  selected  G  for  fused  silica,  and  the  factor  2.8  x  10  3  is  the  ratio  of  permissible  tensile 
stress  to  Young’s  modulus.  Thus 

(A/,) max  “  — 7T  (3 . 1 5  x  10-3)3  x  1.3  X  10716. 

=  7.98  X  10' 5  in  lb. 

Hence  the  maximum  angle  of  twist  per  unit  of  length  is 
7.98  x  10'5  x  32 


fimax 


4.5  x  106  x  7t  x  (3.15  x  \0~3)4 


1.83  rad  105 c 


in. 


in. 


0.72  rad 


cm 


Now  choose  the  angle  of  oblique  incidence  to  be  4°.  Then  for  a  core  which  is  4.5  /u. m  (1.77 
x  1  O'4  in.)  in  diameter 

#.  4sin  9  Mt  x  4  x  sin4° x  7.98  x  10'5 in-lb  (1.77  x  10~4) 

(p  q )  ~~  i  ^  -  i  rrr 


7T  Q 


3.15  x  10 


1-3 


=  102  psi  (703  kPa). 


Hence,  finally,  the  differential  An  due  to  birefringence  is 

An  -  -2.55  x  10'8  x  102  -  -2.60  x  10'6. 

We  see  that  pure  torsion  has  negligible  birefringent  effect  for  a  fiber  carrying  a  single  mode.  Twist, 
however,  introduces  coupled  modes,  even  in  a  single  mode  fiber  (see  the  section  on  coupled  modes). 
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Fiber  Subject  to  Random  Stresses  (Case  8) 

Assume  that  the  stress  held  in  the  fiber  is  randomly  distributed  and  can  be  represented  as  having 
a  transverse  component  of  fluctuation, 

®^rms  ~  &  (p  ~  <?)rms. 

Then  the  effective  rms  difference  between  the  fast  and  slow  components  of  birefringence  is 

From  this  one  obtains  the  rms  phase  difference  between  the  two  components  after  a  length  Ls 
^0  Anrms  =  y/2LLt. 

Our  main  interest  is  to  calculate  Anrms  for  fibers  made  of  fused  silica: 

Anrms  =  -2.55  x  10  885rms  (8s  in  psi) 

1-3.70 5srms  (8s  in  TPa)]. 


We  select  for  calculation  a  value  of  Anrms  which  is  equal  to  that  calculated  for  applied  diametral  forces, 
namely  Anrms  =  —8.5  x  1(T4.  The  corresponding  fluctuation  in  stress  required  to  generate  an  equal  A n 
is 

—  8  S  x  1 0'4 

6srms  *  — — - *  3.3  x  104  psi  (230  MPa) 

rrm  -2.55  x  10“8 


This  is  very  large.  We  conclude  that  average  frozen  fluctuation  stresses  have  negligible  influence  in 
producing  a  An  of  significance  if  they  are  of  the  order  of  1 03 psi  (7  MPa)  or  less.  As  an  extreme  case, 
suppose  An  =  0.1  n  =  0.146.  Then  the  required  frozen  stresses  have  magnitude 

IfiSnJ  =  y55°7^-8  =  5.7  x  1 06 psi  (39GPa). 

Such  stresses  are  far  beyond  the  ultimate  strength  of  the  fiber. 


The  phase  difference  after  length  L  depends  on  Lc. 
An  —  8.5  x  10“4  and  A.  =  5.461  x  1(T7  m  we  have 


Assume  ^/2LL[  =  1  m  =  Leff. 


Then  for 


S<£rms  “ 


27T 


A n rrr,c  L 


rms  ^eff  ' 


2ir  x  8.5  x  10~4  x  } 
5.461  x  1(T7 


=  9.78  x  103rad. 


If  we  desire  50, =  1  rad,  then  8srm5  ~  3.4  psi  (23  kPa)  and  An  —  8.5  x  10  8. 
CONCLUSIONS 


Laser  light  propagating  in  a  fiber  may  experience  a  change  in  the  state  polarization  when  the  fiber 
is  subject  to  mechanical  stress.  The  classical  formula  for  this  effect  is  written  in  terms  of  a  retardation 
in  units  of  wavelengths  (wA0),  or  phase  difference  in  units  of  radians  (A0^)  between  slow  and  fast 
components  of  birefringence: 

m\  =  C(p  -  q)d 


and 


b<t>d  *  k0C(p  -  q)d 
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(symbols  explained  previously  in  text).  An  analysis  of  nine  cases  of  applied  stress  has  been  presented. 
It  is  concluded  that  the  birefringent  effects  are  negligible  in  most  applications  because  there  are  no 
components  of  stress  transverse  to  the  optic  ray.  Two  important  cases  do  show  birefringent  effects: 
fiber  squeezed  by  diametral  forces  and  fiber  subject  to  random  stresses.  Numerical  calculations  predict 
that  for  fibers  of  fused  silica,  at  the  limit  of  allowable  stress, 

A nmdX  —  —8.5  x  10  4  for  the  case  of  applied  diametral  forces 


and 


A <f>d  =  1  rad  if  the  rms  stress  fluctuation  is  3.4  psi  (23  kPa). 

The  case  of  a  fiber  80  /x  m  in  diameter  with  a  core  4.5  /xm  ih  diameter  wound  on  a  mandrel  1  cm  in 
diameter  has  also  been  numerically  estimated  for  fused  silica  at  the  maximum  allowable  stress.  Here 
the  birefringent  effect  is  negligible  if  the  optic  ray  is  parallel  to  the  axis  of  the  fiber.  However  if  the 
optic  ray  is  incident  at  4°  to  the  optic  fiber,  then  it  is  estimated  that 

- 2.5  x  10  4. 
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